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FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC 
EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES 
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Abstract. Over the last ten years, the Einite Element Exterior Calculus (EEEC) has 
been developed as a general framework for linear mixed variational problems, their nu¬ 
merical approximation by mixed methods, and their error analysis. The basic approach 
in EEEC, pioneered by Arnold, Ealk, and Winther in two seminal articles in 2006 and 
2010, interprets these problems in the setting of Hilbert complexes, leading to a more 
general and complete understanding. Over the last five years, the EEEC framework has 
been extended to a broader set of problems. One such extension, due to Holst and Stern 
in 2012, was to problems with variational crimes, allowing for the analysis and nu¬ 
merical approximation of linear and geometric elliptic partial differential equations on 
Riemannian manifolds of arbitrary spatial dimension. Their results substantially gen¬ 
eralize the existing surface finite element approximation theory in several respects. In 
2014, Gillette, Holst, and Zhu extended the EEEC in another direction, namely to para¬ 
bolic and hyperbolic evolution systems by combining the EEEC framework for elliptic 
operators with classical approaches for parabolic and hyperbolic operators, by viewing 
solutions to the evolution problem as lying in Bochner spaces (spaces of Banach-space 
valued parametrized curves). Related work on developing an EEEC theory for parabolic 
evolution problems has also been done independently by Arnold and Chen. In this arti¬ 
cle, we extend the work of Gillette-Holst-Zhu and Arnold-Chen to evolution problems 
on Riemannian manifolds, through the use of framework developed by Holst and Stern 
for analyzing variational crimes. We establish a priori error estimates that reduce to the 
results from earlier work in the flat (non-criminal) setting. Some numerical examples 
are also presented. 
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1. Introduction 

Arnold, Falk, and Winther [2, 3] introduced the Finite Element Exterior Caleulus 
(FEEC) as a general framework for linear mixed variational problems, their numerical 
approximation by mixed methods, and their error analysis. They reeast these problems 
using the ideas and tools of Hilbert complexes, leading to a more complete understand¬ 
ing. Subsequently, Holst and Stern [22] extended the Amold-Falk-Winther framework 
to inelude variational crimes, allowing for the analysis and numerieal approximation of 
linear and geometrie elliptie partial differential equations on Riemannian manifolds of 
arbitrary spatial dimension, generalizing the existing surfaee finite element approxima¬ 
tion theory in several direetions. Gillette, Holst, and Zhu [20] extended the FEEC in 
another direetion, namely to parabolie and hyperbolie evolution systems by eombining 
reeent work on the FEEC for elliptie problems with a elassieal approaeh of Thomee [33] 
to solving evolution problems using semi-diserete finite element methods, by viewing so¬ 
lutions to the evolution problem as lying in Boehner spaees (spaees of Banaeh-spaee val¬ 
ued parametrized eurves). Arnold and Chen [1] independently developed related work, 
for generalized Hodge Eaplaeian parabolie problems for differential forms of arbitrary 
degree, and Holst, Mihalik, and Szypowski [21] eonsider similar work in adaptive finite 
element methods. In this artiele, we aim to combine the approaches of the above articles. 
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extending the work of Gillette, Holst, and Zhu [20] and Arnold and Chen [1] to para- 
bolie evolution problems on Riemannian manifolds by using the framework of Holst and 
Stem [22]. 

1.1. The Hodge heat equation and its mixed form. We now introduee our problem 
by some eonerete motivation. We eonsider an evolution equation for differential forms 
on a manifold. Then we rephrase it as a mixed problem as an intermediate step toward 
semidiseretization using mixed finite element methods. We then see how this allows us 
to leverage existing a priori error estimates for parabolie problems, and see how it fits in 
the framework of Hilbert eomplexes. 

Let M be a eompact oriented Riemannian n-manifold embedded in The Hodge 
heat equation is to find time-dependent fc-form m : M x [0, T] —)■ A^(M) sueh that 

Ut — Au = Ut + {Sd + d6)u = f in M, for t > 0 

u{-, 0) = g in M. ^ ' 

where g is an initial fc-form, and /, a possibly time-dependent fc-form, is a source term. 
Note that no boundary conditions are needed for manifolds without boundary. This is 
the problem studied by Arnold and Chen [1], and in the case k = n, one of the problems 
studied by Gillette, Holst, and Zhu [20], building upon work in special cases for domains 
in and by Johnson and Thomee [24, 33] 

For the stability of the numerical approximations with the methods of [23] and [3], we 
recast the problem in mixed form, converting the problem into a system of differential 
equations. Motivating the problem by setting a = 5u (recall that for the Drrichlet prob¬ 
lem and k = n, 5 here corresponds to the gradient in Euclidean space, and is the adjoint 
d, corresponding to the negative divergence), and taking the adjoint, we have 

{a,uj) — {u,du) =0, V ce G f > 0, 

{ut,(p) + {da,ip) + {du,dip) ={f,ip), \/ p E HVL^{M) f > 0. (1.2) 

m(0) = g. 

Unlike the elliptic case, we do not have to explicitly account for harmonic forms in the 
formulation of the equations themselves, but they will definitely play a critical role in our 
analysis and bring new results not apparent in the k = n case. 

1.2. Semidiscretization of the equation. In order to analyze the numerical approxima¬ 
tion, we semidiscretize our problem in space. In our case, we shall assume, following 
[22], that we have a family of approximating surfaces Mh to the hypersurface M, given 
as the zero level set of some signed distance function, all contained in a tubular neigh¬ 
borhood U of M, and a projection a : Mh —)■ M along the surface normal (of M). The 
surfaces may be a triangulations, i.e., piecewise linear (studied by Dziuk and Demlow in 
[16, 14]), or piecewise polynomial (obtained by Lagrange interpolation over a triangu¬ 
lation of the projection a, as later studied by Demlow in [13]). We pull forms on Mh to 
M back via the inverse of the normal projection, which furnishes injective morphisms 

HVL^{M) as required by the theory in [22], which we shall review in Section 
2 below. Finally, we need a family of linear projections 11^ : HVL^{M) —>■ A^ such 
that flh° ih = id which allow us to interpolate given data into the chosen finite element 
spaces—this is necessary because some of the more obvious, natural seeming choices 
of operators, such as i\, can be difficult to compute (nevertheless, will still be useful 
theoretically). 
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We now can formulate the semidiscrete problem: we seek a solution {ah, Uh) € Hh x 
Sh C X HVt^ such that 

{ah,ujh)h-{uh,duh)h =0, WuhEHh, t>0 

{uh,t,^h)h +{dah,(Ph)h +{duh,diph)h = {Ilhf,iph)h, (Ph ^ Su f>0 (1.3) 

Uh{0) = gh- 

We shall describe how to define gh € Sh shortly; it is to be some suitable interpolation 
of g. As Sh and Hh are finite-dimensional spaces, we can reduce this to a system of 
ODEs in Euclidean space by choosing bases {^jJi) for Sh and (0^) for Hh', expanding the 
unknowns ah = and Uh = U^{^)4>k', substituting these basis functions 

as test functions to form matrices Am = {(j)k, (j)e), Bik = {dijji, (j)k), and Dij = {'ipi, ipj)', 
and finally forming the vectors for the load data F defined by F*. = (F, 0*.), and initial 
condition G defined hy gh = We thus arrive at the matrix equations for the 

unknown, time-dependent coefficient vectors S and U: 

DE-B^U = 0 , 

AUt + BE + KU = F, forf>0 
(7(0) = G. 

The matrices A and D are positive definite, hence invertible. Substituting S = D~^B^U, 
we have the system of ODEs 

AUt + {BD-^B^ + K)U = F, for t > 0, f/(0) = G, 


which has a unique solution by the usual ODE theory. Eor purposes of actually numeri¬ 
cally integrating the ODE, namely, discretizing fully in space and time, it is better not to 
use the above formulation, because it can lead to dense matrices. Computationally, this 
is due to the explicit presence of an inverse, D~^, not directly multiplying the variable; 
conceptually, this is actually a statement about the discrete adjoint to the codifferential 
dj^ generally having global support even if the finite element functions are only locally 
supported [1]. Instead, we differentiate the first equation with respect to time, getting 
DTit — B^Ut = 0, which leads to the block system 


d 

dt 


D 

0 


-B^ 

A 


0 

-B 


0 

-K 


+ 


(1.4) 


which is still well-defined ODE for S and [/, as the invertible matrices A and B appear 
on the diagonal. This differentiated equation also plays a role in the showing that the 
continuous problem is well-posed. 

These equations differ from those studied by Gillette, Holst, and Zhu [20], Arnold 
and Chen [1], and Thomee [33] by the choice of finite element spaces—here we are 
assuming them to be in some Sobolev space of differential forms on manifolds (or in 
a triangulated mesh in a tubular neighborhood) rather than subsets of Euclidean space. 
This suggests that we should try to gather these commonalities, examine what happens 
in abstract Hilbert complexes, and see how general a form of error estimate we can get 
this way. 


1.3. Error analysis. The general idea of the method of Thomee [33] is to compare the 
semidiscrete solution to an elliptic projection of the data, a method first explored by 
Wheeler [35]. If we assume that we already have a solution u to the continuous problem, 
then for each fixed time t, u{t) can be considered as trivially solving an elliptic equation 
with data —Au{t). Thus, using the methods developed in [3], we consider the discrete 
solution Uh for u in this elliptic problem (namely, applying the discrete solution operator 
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T/j to —Au{t)). This may be compared to the true solution (at each fixed time) using 
the error estimates in [3] . What remains is to compare the semidiscrete solution Uh (as 
defined by the ODEs (1.3) above) to the elliptic projection, so that we have the full error 
estimate by the triangle inequality. Thomee derives the following estimates, for finite 
elements in the plane (n = 2) of top-degree forms (k = 2, there represented by a scalar 
proxy), for gh the elliptic projection of the initial condition g and t > 0: 

\\Uh{t) -u{t)\\L2 ^\\u{t)\\H2 + ||Mi(s)||H2ds^ , (1.5) 


o-hit) 


CT(f)||L2 < ch^ 


u{t)\\H3 + 



( 1 . 6 ) 


Gillette, Holst, and Zhu [20], and Arnold and Chen [1] generalize these estimates and 
represent them in terms of Bochner norms. These estimates describe the accumulation 
of error up to fixed time value t, assuming, of course, that the spaces finite elements 
are sufficiently regular to allow those estimates. The key equation that makes these 
estimates possible are Thomee’s error evolution equations: defining p = \\uh{t) — M(f)||, 
6 = \\uh{t) - Uh{t)\\, and e = ||a/j(f) - dh(f)||, we have 


{e,Uh) + {6,divujh) = 0 . 


These are used to derive certain differential inequalities and make Gronwall-type esti¬ 
mates. In this chapter, we examine the above error equations and place them in a more 
abstract framework. We use Bochner spaces (also used by [20]) to describe time evo¬ 
lution in Hilbert complexes, building on their successful use in elliptic problems. We 
investigate Thomee’s method in this framework to gain further insight into how finite 
element error estimates evolve in time. 


1.4. Summary of the paper. The remainder of this paper is structured as follows. In 
Section 2, we review the finite element exterior calculus (FEEC) and the variational 
crimes framework of Holst and Stern [22]. We prove some extensions in order to ac¬ 
count for problems with prescribed harmonic forms; this is what allows the elliptic pro¬ 
jection to work in the case where harmonic forms are present. In Section 3, we formulate 
abstract parabolic problems in Bochner spaces and extend some standard results on the 
existence and uniqueness of strong solutions, and describe how this problem fits into 
that framework. In Section 4, we extend the a priori error estimates for Galerkin mixed 
finite element methods to parabolic problems on Hilbert complexes. Then, we relate 
the resuls to the problem on manifolds. The main abstract result is Theorem 4.2, which 
uses the previous results from the FEEC framework with variational crimes, in order to 
understand how those error terms evolve with time. We then specialize, in Section 5 to 
parabolic equations on Riemannian manifolds, our original motivating example, and see 
how this generalizes the error estimates of Thomee [33], Gillette, Holst, and Zhu [20], 
and Holst and Stem [22]. In Section 6, we present a numerical experiment comparing 
the methods based on this mixed form to more straightforward implementations in the 
scalar heat equation case. 


2. Finite Element Exterior Calculus 

We review here the relevant results from the finite element exterior calculus (FEEC) 
that we will need for this paper. FEEC was introduced in Arnold, Falk and Winther [2, 3] 
as a framework for deriving error estimates and formulating stable numerical methods for 
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a large class of elliptic pde. One of the central ideas which helped unify many of these 
distinct methods into a structured framework has been the idea of Hilbert complexes [8], 
which abstracts the essential features of the cochain complexes commonly found in exte¬ 
rior calculus and places them in a context where modern methods of functional analysis 
may be applied. This assists in formulating and solving boundary value problems, in 
direct analogy to how Sobolev spaces have helped provide a framework for solving such 
problems for functions. Arnold, Falk, and Winther [3] place numerical methods into this 
framework by choosing certain finite-dimensional subspaces satisfying certain compat¬ 
ibility and approximation properties. Holst and Stem [22] extended this framework by 
considering the case in which there is an injective morphism from a finite-dimensional 
complex to the complex of interest, without it necessarily being inclusion. This allows 
the treatment of geometric variational crimes [6, 7], where an approximating manifold 
(on which it may be far easier to choose finite element spaces) no longer coincides with 
the actual manifold on which we seek our solution. We review the theory as detailed in 
[22] and refer the reader there for details. 

2.1. Hilbert Complexes. As stated before, the essential details of differential com¬ 
plexes, such as the de Rham complex, are nicely captured in the notion of Hilbert com¬ 
plexes. This enables us to see clearly where many elements of boundary value problems 
come from, in particular, the Laplacian, Hodge decomposition theorem, and Poincare 
inequality. In addition, it allows us to see how to carry these notions over to numerical 
approximations. 

Definition 2.1 (Hilbert complexes). We define a Hilbert complex (IF, d) to be sequence 
of Hilbert spaces with possibly unbounded linear maps : C*’ C IF^ —)■ C 

IF*^’''^, such that each d^ has closed graph, densely defined, and satisfies the cochain 
property d^ o d^~^ = 0 (this is often abbreviated d? = 0,' we often omit the superscripts 
when the context is clear). We call each the domain of d^. We will often refer to 
elements of such Hilbert spaces as ‘forms,” being motivated by the canonical example 
of the de Rham complex. The Hilbert complex is called a closed complex if each image 
space 53^ = (called the k-coboundaries is closed in W^, and a bounded 

complex if each df is in fact a bounded linear map. The most common arrangement in 
which one finds a bounded complex is by taking the sequence of domains V^, the same 
maps df, but now with the graph inner product 

{v,w)v = {v,w) -I {d^v,d^w). 

for all v,w E V^. Unsubscripted inner products and norms will always be assumed to 
be the ones associated to W^. 

Definition 2.2 (Cocycles, Coboundaries, and Cohomology). The kernel of the map d^ 
in will be called 3^, the k-cocycles and, as before, we have = d^~^V^~^. Since 
dk Q df-i _ have C 3^, so we have the k-cohomology 3^/®^- The harmonic 
space is the orthogonal complement of^^ in 3^. This means, in general, we have 
an orthogonal decomposition and we have that is isomorphic to 

3^/®^, the reduced cohomology, which of course corresponds to the usual cohomology 
for closed complexes. 

Definition 2.3 (Dual complexes and adjoints). For a Hilbert complex (IF, d), we can 
form the dual complex (IF*, d*) which consists of spaces IF^ = IF^, maps df : Vf C 
IF^ —)■ C IF^_^ such that d\j^^ = (d^)*, the adjoint operator, that is: 

{dl^iv,w) = {v,dl"w). 
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The operators d* decrease degree, so this is a chain complex, rather than a cochain 
complex; the analogous concepts to cocycles and coboundaries extend to this case and 
we write 3^ and ^Ifor them. 

Definition 2.4 (Morphisms of Hilbert complexes). Let {W, d) and [W, d') be two Hilbert 
complexes, f : W ^ W is called a morphism of Hilbert complexes if we have a 
sequence of bounded linear maps —)■ W'^ such that d"^ o /^ = o d’^ (they 

commute with the differentials). 

With the above, we can show the following weak Hodge decomposition: 

Theorem 2.5 (Hodge Decomposition Theorem). Let (W, d) be a Hilbert complex with 
domain complex (V, d). Then we have the W- and V-orthogonal decompositions 

= ^ © © 3^^w (2.1) 

. (2.2) 

where n V^. 

Of course, if 03^ is closed, then the extra closure is unnecessary, and we omit the term 
“weak”. We shall simply write for , which is will be the most useful orthogonal 
complement for our purposes. The orthogonal projections Pu for a subspace Lf will be 
in the Vk-inner product unless otherwise stated (although again, due to the two inner 
products coinciding on 3^ and its subspaces, they may be the same). We note that by the 
abstract properties of adjoints ([3, §3.1.2]), and = 3^- Also very 

useful is that the V- and PT-norms agree on 3 and hence on and i^. 

The following inequality is an important result crucial to the stability of our solutions 
to the boundary value problems as well as the numerical approximations: 

Theorem 2.6 (Abstract Poincare Inequality). If(V, d) is a closed, bounded Hilbert com¬ 
plex, then there exists a constant cp > 0 such that for all v G 

ll^'llv < cp\\d^v\\v- 

In the case that (V, d) is the domain complex associated to a closed Hilbert complex 
(W,d), {V,d) is again closed, and the additional graph inner product term vanishes: 

= ||(i^n||. We now introduce the abstract version of the Hodge Laplacian and the 
associated problem. 

Definition 2.7 (Abstract Hodge Laplacian problems). We consider the operator L = 
dd* + d*dona Hilbert complex (W, d), called the abstract Hodge Laplacian. Its domain 
is Dl = {m G n Vf : du G d*u G and the Hodge Laplacian problem 

is to seek u ®Vk, given f G W^, such that 

{du,dv) + {d*u,d*v) = {f,v) (2.3) 

for all u G n Vf. This is simply the weak form of the Laplacian and any m G fl Vf 

satisfying the above is called a weak solution. Owing to difficulties in the approximation 
theory for such a problem (it is difficult to construct finite elements for the space fl I 4 *), 

Arnold, Falk, and Winther [3] formulated the mixed abstract Hodge Laplacian problem 
by defining auxiliary variables a = d*u and p = Psjf, the orthogonal projection of 
f into the harmonic space, and considering a system of equations, to seek (a,u,p) G 
X X such that 

{a,T) — (u^dr) =0 Vr G 
{da, v) + {du, dv) + {p, v) = (/, v) \/v G 
{u, g) = 0 Vg G 


(2.4) 
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The first equation is the weak form of a = d*u, the seeond is the main equation 
(2.3) modified to aeeount for a harmonie term so that a solution exists, and the third 
enforces uniqueness by requiring perpendicularity to the harmonic space. With these 
modifications, the problem is well-posed by considering the bilinear form (writing := 
X X S)^) S : X —)■ M defined by 

B{a,u,p;T,v,q) := (cr,r) — {dT,u) + {da,v) + {du,dv) + {p,v) — {u,q). (2.5) 

and linear functional F G (j£^)* given by F(r, v, q) = (/, v). The form B is not coercive, 
but rather, for a closed Hilbert complex, satisfies an inf-sup condition [3, 4]: there exists 
7 > 0 (the stability constant) such that 

. f B{a,u,p-,T,v,q) 

mf sup 777 -TT—f 77 -^ =: 7 > 0 . 

{a,u,p)j^0 (r,v,q)^0 || P) ||t|| (T V, Q) llx 

where we have defined a standard norm on products: ||(cr, u,p)||x := ||<^||y + ||'?^||y + ||p||- 
This is sufficient to guarantee the well-posedness [4]. To summarize: 

Theorem 2.8 (Arnold, Falk, and Winther [3], Theorem 3.1). The mixed variational prob¬ 
lem (2.4) on a closed Hilbert complex (kF, d) with domain (V, d) is well-posed: the bi¬ 
linear form B satisfies the inf-sup condition with constant 7 , so for any F G (;T^)*, 
there exists a unique solution (a, u, p) to (2.4), i.e., B{a, u, p; r, v, q) = F(r, u, q) fo all 
(r, u, q) G and moreover, 

Wia,u,p)\\x < 7"ii^|k*- 

The stability constant 7“^ depends only on the Poincare constant. 

Note that the general theory (e.g., [4, 18]) guarantees a unique solution exists for any 
bounded linear functional F G (X^)*, which in this case with product spaces, means that 
the problem is still well-posed when there are other nonzero linear functionals on the 
RHS of (2.4) besides {/, v). We shall need this result for parabolic problems, where we 
assume u has a harmonic part (Pf^u 7 ^ 0 ). 

2.2. Approximation of Hilbert Complexes. We now approximate solutions to the ab¬ 
stract mixed Hodge Laplacian problem. To do so, Arnold, Falk, and Winther [3] in¬ 
troduce finite-dimensional subspaces 14 C 1 / of the domain complex, such that the 
inclusion 4 : 14 C is a morphism, i.e. dV^ C With the weak form (2.4), we 

formulate the Galerkin method by restricting to the subspaces: 

{ah, t) - {uh, dr) =0 Vr G 

{dah,v) + {duh,dv) + {ph,v) =(/,u) Vu G (2.6) 

{Uh,q) =0 VgG.^t 

We abbreviate by setting X^ := x x We must also assume the existence 
of bounded, surjective, and idempotent (projection) morphisms iih '■ V —)■ 14- It is 
generally not the orthogonal projection, as that fails to commute with the differentials. 
As a projection, it gives the following quasi-optimality result: 

||m — 7thu\\v = inf IK-I — '^h){u — n)||\/ < ||/ — 7ih\\ inf 14 — u||\/. 

veVh veVh 

The problem (2.6) is then well-posed, with a Poincare constant given by cp|| 7 r^||, where 
cp is the Poincare constant for the continuous problem. This guarantees all the previous 
abstract results apply to this case. With this, we have the following error estimate: 
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Theorem 2.9 (Arnold, Falk, and Winther [3], Theorem 3.9). Let (14, d) be a family of 
subcomplexes of the domain (V, d) of a closed Hilbert complex, parametrized by h and 
admitting uniformly V-bounded cochain projections Tih, cind let {a,u,p) G be the 
solution of the continuous problem and {ah, Uh,Ph) G ^h corresponding discrete 

solution. Then the following error estimate holds: 

\\{a - ah, u - Uh,p - Ph)\\x = \W - CThWv + \\u - Uh\\v + \\p - Ph\\ 

< C{ inf ||(T — t\\v+ inf ||m — f||\/ + inf \\p — q\\v + b inf \\P^u — nHy) 

’'S14 9^14 ’'64? 

(2.7) 

with p = p^ = suprg^fe II (7 — TT^) r ||, the operator norm of I — restricted to 

Corollary 2.10. If the I 4 approximate V, that is, for all u E V, inf.ygy^ ||m — n||y —)■ 0 
as h ^ 0, we have convergence of the approximations. 

In general, the harmonic spaces and do not coincide, but they are isomorphic 
under many circumstances we shall consider (namely, the spaces are isomorphic if for all 
harmonic forms q G the error \\q — TThq\\ is at most the norm ||g|| itself [3, Theorem 
3.4], and it always holds for the de Rham complex). For a quantitative estimate relating 
the two different kinds of harmonic forms, we have the following 

Theorem 2.11 ([3] .Theorem 3.5). Let (V, d) be a bounded, closed Hilbert complex, 
( 14 , d) a Hilbert subcomplex, and iih a bounded cochain projection. Then 

ll(/ - -R)j9lk < ll(/ - e ( 2 . 8 ) 

ll(/-A,)<)llv< ll(/-4)T)9lk,Vgefit (2.9) 

2.3. Removing the Suhcomplex Assumption: Variational Crimes. For geometric 
problems, it is essential to remove the requirement that the approximating complex I 4 
actually be subspaces of V. This is motivated by the example of approximating planar 
domains with curved boundaries by piecewise-linear approximations, resulting in finite 
element spaces that lie in a different function space [6]. Holst and Stem [22] extend the 
Arnold, Falk, Winther [3] framework by supposing that ih ■ Vh ^ V is an injective 
morphism which is not necessarily inclusion; they also require projection morphisms 
vT/i : V —)■ 14 with the property tt/j o 4 = id, which replaces the idempotency re¬ 
quirement of the preceding case. To summarize, given {W, d) a Hilbert complex with 
domain (V, d), (H4, dh) another complex (whose inner product we denote (-, ■)h) with 
domain (I 4 , dh), injective morphisms 4 : W 4 W, and finally, projection morphisms 
Tih'.V —)■ 14 . We then have the following generalized Galerkin problem: 

{cTh, Th)h - {uh, dhThjh =0 Vt/, G 14 ^"^ 

{dhah,Vh)h + {dhUh,dhVh)h + {Ph,Vh)h = {fh,Vh)h G V/j (2.10) 

{uh,qh)h =0 WqhEft’j, 

where fh is some interpolation of the given data / into the space Wh (we will discuss 
various choices of this operator later). This gives us a bilinear form 

Bh {^h, '^h, Ph, Th, Vh, qh) ■ {^h, '^h}h {^h, dhTh) h 

+ {dhah,Vh)h + {dhUh,dhVh)h + {Ph,Vh)h - {uh,qh)h- (2.11) 

This problem is well-posed, which again follows from the abstract theory as long as the 
complex is closed, and there is a corresponding Poincare inequality: 
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Theorem 2.12 (Holst and Stem [22], Theorem 3.5 and Corollary 3.6). Let {V,d) and 
(14, dh) be bounded closed Hilbert complexes, with morphisms ih ■ Vh ^ V and Tih : 
1/ —)■ 14 such that vT/i o 4 = id. Then for all Vh € we have 

\\vh\\vh < cp llvr^ll \\dhVh\\vh, 

where cp is the Poincare constant corresponding to the continuous problem. If {V,d) 
and (14, dfi) are the domain complexes of closed complexes (14, d) and {Wh, df), then 
Ild/iL/illvh i.s simply (since it is the graph norm and df = D). 

In other words, the norm of the injeetive morphisms 4 also eontributes to the stability 
eonstant for this diserete problem. Analysis of this method results in two additional 
error terms (along with now having to explicitly reference the injective morphisms 4 
which may no longer be inclusions), due to the inner products in the space 14 no longer 
necessarily being the restriction of that in V: the need to approximate the data /, and the 
failure of the morphisms 4 to be unitary. 

Theorem 2.13 (Holst and Stern [22], Corollary 3.11). Let (V, d) be the domain complex 
of a closed Hilbert complex (14, d), and (Vh, dh) the domain complex of (Wh, dh) with 
morphisms ih ■ Wh —)■ 14 and 7ih '■ V ^ Vh as above. Then if we have solutions (a, u,p) 
and ((Th,Uh,Ph) to (2.4) and (2.10) respectively, the following error estimate holds: 

||o- - 4o-h||y + \\u - ihUhWv + lb - ihPhW 

<C( inf Ik —r||y+ inf |b —r;||y+ inf \\p — q\\v + h inf \\P<su — v\\v 

T&hV^-^ veihV^ q&ihV^ 

+ IIA-*^/IU + ||/-411 ll/ll), (2.12) 

where Jh = 44, and p = p,^= sup || (f — ih'^h) ^||- 

The extra terms (in the third line of the inequality) are analogous the terms described 
in the Strang lemmas [ 6 , §III.l]. The main idea of the proof of Theorem 2.13 (which we 
will recall in more detail below, because we will need to prove a generalization of it as 
part of our main results) is to form an intermediate complex by pulling the inner products 
in the complex (14, d) back to (Wh, dh) back by 4, construct a solution to the problem 
there, and compare that solution with the solution we want. This modified inner product 
does not coincide with the given one on Wh precisely when 4 is not unitary: 

{v,w)i*^w = {ihV,ihw)h = {ilihV,w)h = {JhV,w)h. 

Unitarity is simply the condition Jh = I. The complex Wh with the modified inner 
product now may be identified with a true subcomplex of 14, for which the theory of 
[3] directly applies, yielding a solution (a'j^,u'j^,p'f} G 4^^“^ x V!^ x where 
is the discrete harmonic space associated to the space with the modified inner prod¬ 
uct. This generally does not coincide with the discrete harmonic space since the 
discrete codifferential d*f in that case is defined to be the adjoint with respect to the 
modified inner product, yielding a different Hodge decomposition. The estimate of 
1144 “ *^14 + Il 44 ~ '^llv + ll 4 Ph — pII then proceeds directly from the preceding 
theory for subcomplexes (2.7). The variational crimes, on the other hand, arise from 
comparing the solution (<Jh,Uh,Ph) with (a'j^,u'j^,p'f). Finally, the error estimate (2.12) 
proceeds by the triangle inequality (and the boundedness of the morphisms 4 )- 
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2.4. Elliptic Error Estimates for a Nonzero Harmonic Part. Our objective in the 
remainder of this section is to prove one of our main results, a generalization of Theorem 
2.13 which allows the possibility of the solution u having a nonzero harmonic part w. 
We first need a few lemmas. 

Lemma 2.14. Theorem 2.9 continues to apply when we have {u, q) = (w, q) for all 
q G where w E is prescribed (i.e., = w, which may generally not be zero). 

Proof. We closely follow the proof, in [3], of Theorem 2.9 above, noting where the mod¬ 
ifications must occur. Let B be the bounded bilinear form (2.5); then {a,u,p) satisfies, 

for all {Th,Vh,qh) e 

B{a,u,p;Th,Vh,qh) = {f,Vh) - {u,qh). 

We (/-orthogonally project {a,u,p) in each factor to {T,v,q) G X^. Then for any 

{Th,Vh,qh) G X^, 

B{ah -T,Uh- V, ph - q] Th, Vh, qh) 

= B{a - T,u - v,p - q]Th,Vh,qh) + {u,qh) - {w,qh) 

= B{a -T,u-v,p- q;Th,Vh,qh) + {Psj^{u - w), qh) 

< C (||cr - t\\v + \\u - nlly -f Up - g|| 4- \\Psj^{u - m;)||) {\\Th\\v + ||^^/i||y + ||g/ii|). 

(2.13) 

Noticing that the factor ph — q in the bilinear form above is in the original domain 
we can now choose the appropriate (r/i, Vh,qh) that verifies inf-sup condition of B: 

B{(Th -T,Uh- V, Ph - q; Th, Vh, qh) 

> yiWah - t\\v + \\uh - v\\v + \\ph - glDdlThHu + \\vh\\v + lkh||)- 

Comparing this to (2.13) above, we may cancel the common factor, and divide by 7 to 
arrive at 


\Wh -t\\v + \\Uh -v\\v + \\ph - g|| 

< (||cr - t\\v + ||m - v\\v + \\p - q\\ + \\Pf}^,{u - «;)||) . (2.14) 

This differs (aside from the notation) from [3] in that we now have, rather than P^j^u, 
instead P^^{u — w), with the harmonic part subtracted off. Removing the harmonic part 
allows us to continue as in [3]: the Hodge decomposition u — w = u — Ps^u consists only 
of coboundary and perpendicular terms + u±_ G 93^ © With contained in 
3^, it follows Ps)^u± = 0, and = 0. Also, (/ — Th)u<s is perpendicular to 

Therefore, for dXXq E 


- Psju), q) = {Psj^u^, q) = q) 

= (w® — ThUfs, q) = (w® — ThU<s, (/ — Pffq) ■ 

Now, setting 

we have 


WP^ihiP - P^'^)\\ = {Pfyni^ - Pf:u),(l) = (m® - ThU<B, (/ - Pf,)q) 

< \\u<B - ThU<B\\ ||(/ - Pr-))q\\ < C\\{I - Psj)q\\ inf ||m® - 
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Finally, by the second estimate of Theorem 2.11 above, we can bound ||(/ — Pfj)q\\ by 
11(1 - 7rh)Psjq\\, giving us 

II (I - PfiM < II (^ - ^h)Pf}q\\ < sup^ II (/ - nh)r\\ \\Pf,q\\ < /i. 

From the triangle inequality, we derive the estimate 

llo- - ahWv + \\u - UhWv + lb - PhW 

< lb - rib + lb - v\\v + lb - qW + lb - <^h\\v + \\uh - v\\v + lb -phll 

< (1 + C'7"b ( Ik ~'r||y + lb ~+ lb “ ?ll +/^ \\P^u-v\\v 

V 

Using best approximation property of orthogonal projections, we can express the remain¬ 
ing terms with the infima, and this gives the result. □ 

We also need a technical lemma which enables us to identify the orthogonal projection 
onto the identified subcomplex ih^h order to be able to make additional estimates 
of the variational crimes in terms of the operator norms ||J — J/j||. It is the infinite¬ 
dimensional analogue of taking a Moore-Penrose pseudoinverse [32, §3.3] for infinite¬ 
dimensional spaces: 

Lemma 2.15. Let ih ■ Wh W be an injective map of Hilbert spaces, and J = ifih- 
Then Jh is invertible, and is the Moore-Penrose pseudoinverse ofih, i-c. it maps 

ihWh isometrically back to Wh with the modified inner product. 

We write k foi" Jh^ih- 

Proof. The invertibility of Jh follows directly from the injectivity of ih, which makes 
{Jh-, ■)h a positive-definite form. Now, {Jf^ifijih = Jf^Jh = idw,,, which shows that 
it is in fact a left inverse, as required for pseudoinverses. To show the orthogonality, 
minimizing \\\ihW — h\\^ for any h yields, by the completeness of Wh, the solution 
w = Jfi^i*yjb, showing that it is a least squares solution, therefore the Moore-Penrose 
pseudoinverse. □ 

We are now ready to prove our main elliptic error estimate, an extension of Theorem 
2.13. 

Theorem 2.16 (Extension of elliptic error estimates to allow for a harmonic part). Con¬ 
sider the problems (2.4) and (2.10) but instead with now with prescribed, possibly nonzero 
harmonic part w: Given f G W^ and w G we seek {a, u, p) G such that 

{a, t) — {u, dr) =0 Vr G 

{da,v)-\-{du,dv)-\-{p,v) = {f,v) Vl G (2.15) 

{u, q) = {w, q) Vg G 

The solution to this problem exists and is unique, with w indeed equal to P^u, and is 
bounded c(||/|| -f |k||), with c depending only on the Poincare constant. Now, con¬ 
sider the discrete problem, with fh,Wh G 

(cTh, Th)h - {Uh, dhTh)h =0 Vr/i G vlj~^ 

{dhCTh, Vh)h + {dhUh, dhVh)h + {Ph, Vh)h = {fh, Vh)h 'Jvh G Vjf 

{uh,qh)h = {wh,qh)h Wqhefil. 



( 2 . 16 ) 
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This problem is also well-posed, with the modified Poincare constant in Theorem 2.12. 
Then we have the following generalization of the error estimate (2.12) above: 

||(T - ih(Th\\v + \\u - ihUhWv + \\p - ihPhW 

<(7 inf ||(T —r||y+ inf ||m —t>||y+ inf \\p — q\\v + P inf \\P<su — v\\v 

\r£ihV/;;~^ q&ihV^ vei^V^ 

+ inf ilw-^llv + \\fh-i*hf\\h + \\wh-ilw\\h + ||/- A|| (||/|| + HwH) | , 

(2.17) 

where, as before, Jh = and /i = = sup || (7 — ih'^h) ’"||- 

We see that three new error terms arise from the approximation of the harmonie part, 
one being the data interpolation error (but measured in the 14 -norm, whieh partially 
eaptures how d fails to eommute with and how Wh may not neeessarily be a diserete 
harmonie form), another best approximation term, and finally another term from the non- 
unitarity. The relation of fh to / and Wh to w need not be further speeified, beeause the 
theorem direetly expresses sueh a dependenee in terms of their relation to 4/ and ilw; 
it has been isolated as a separate issue. However as mentioned in the introduetion, and 
following [22], we often take fh = H/i/, where lih is some family of linear interpolation 
operators with II/j o 4 = id. Another seemingly obvious ehoiee is 4 itself (thus making 
those eorresponding error terms zero), but as mentioned in [22], this ean be diffieult to 
eompute, so we do not restriet ourselves to this ease. Various ehoiees of interpolation 
will be erueial in deeiding whieh estimates to make in the parabolie problem. We split 
the proof of this theorem into two parts, the first of whieh derives the quantities on the 
seeond line of (2.17), and the seeond part, we derive the quantities on the third line of 
(2.17). Generally, we follow the pattern of proof in [3, Theorem 3.9] and [22, Theorem 
3.10], noting the neeessary modifieations, as well as a similar teehnique given for the 
improved error estimates by Arnold and Chen [1]. 

First part of the proof of Theorem 2.16. First, following Holst and Stern [22] as above, 
we eonstruet the eomplex Wh but with the modified inner produet {Jh'r) (the assoeiated 
domain eomplex 14 remains the same). This gives us a diserete Hodge deeomposition 
with another type of orthogonality and eorresponding diserete harmonie forms and or¬ 
thogonal eomplement (due to a different adjoint df): 

© 3 ^' 

(generally, primed objeets will represent the eorresponding objeets defined with the mod¬ 
ified inner produet; the diserete eoboundaries are in faet the same as before, beeause d 
and dh do not depend on the ehoiee of inner produet). The main eomplieations arise in 
having to keeping eareful traek of the different harmonie forms involved, beeause their 
nonequivalenee and possible non-preservation by the operators eontribute direetly to the 
error. We then define, as in [22], the intermediate solution (a]^, p'f) G x Vjf x 
(whieh we abbreviate as 

{Jhfj'h, n)h - {Jhu'h, dhTh)h =0 Vr^ e 

{Jhdh(j'h,Vh)h + {Jhdhu'h,dhVh)h + {JhPhWh)h = {i*hf,Vh)h G (2.18) 

{Jhu'hW'h)h ={ilw,q'f)h 


14 


M. HOLST AND C. TIEE 


and the corresponding bilinear form x —)■ M given by 

Bki^h^Uh^Ph^^h.yh.Qh) — {Jh(^h^rh)h-{JhUh,dhTH)h 

{Jndh^h") '^h}h {Jhdh'^fi') dh'^h^h ^JhPhi h Qh)h' (2-19) 

This satisfies the inf-sup condition with Poincare constant cp||7r/i||. Note that we will 
need to extend all the bilinear forms B^, and in the last factor to all of V/^ in or¬ 
der to compare the two, since they are initially only defined on the respective, differing 
harmonic form spaces. This is not a problem so long as we remember to invoke the inf- 
sup condition only when using the non-extended versions. The idea is, again, to use the 
triangle inequality: 

||cT - ih(^h\\v + Ik - ihThWv + lb - ihPhW < (2.20) 

lb - ih(rh\\v + Ik - ihThWv + lb - ihPhW (2-21) 

+ \\ih{(Th - (^h)\\v + \\ih{Th - Th)\\v + \\ih{p'h-Ph)\\- (2.22) 

These quantities can be estimated using only geometric properties of the domain; we 
have no need to actually explicitly compute (akrikK)- To estimate the term (2.21) 
(which we shall refer to as the PDE approximation term, whereas (2.22) will be called 
variational crimes), we recall that ih is an isometry of Wh with the modified inner product 
to its image, which is a subcomplex. 

Thus, Lemma 2.14 above applies, with the approximation (4crk ihP'h) on identi¬ 
fied subcomplex ih^h- This gives us the terms on the second line of (2.17). □ 

To finish our main proof, we need to consider the variational crimes (2.22). Since the 
maps ih are bounded, and we eventually absorb their norms into the constant C above, it 
suffices to consider \\ah — o-'hWvh + lb/* ~ '^klt4 + \\Ph ~ p'hWh^ which we shall state as a 
separate theorem. 

Theorem 2.17. Let {<Jh,Uh,Ph) & X’l be a solution to (2.16), u'f^, p'f^) G a solu¬ 

tion to (2.18), and w = P^u, the prescribed harmonic part of the continuous problem. 
Then 


^h-(^'h\\v^ + \\Uh-u'h\\vn + \\Ph-p'h\\h 
< C{\\fh - ilfWh + \\wh - ilw\\v^ + Ik - A||(||/|| + Ibll) + inf 


w^-^lk), 

(2.23) 


i.e., they are bounded by the terms on the third line in (2.17). 

Proof of Theorem 2.17 and second part of the proof of Theorem 2.16. We follow the proof 
of Holst and Stern [22, Theorem 3.10] and note the modifications. Let (r, u,g) and 
(r/i, Vh, Qh) £ Xk Consider the bilinear form Bh, (2.11) above, and write 


Bh{(rh -T,Uh- v,ph - q] Th, Vh,Wh) = Bhffh - a'h, Uh - u'h,Ph - p'h] Vh, Qh) 

+ Bh{(T'h - L u'h - v,p'i^ - g; t/,, vh, Qh)- 
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We then have, recalling the modified bilinear form B'j^, (2.19) above, and extending it in 
the last factors to all of , 

BhWh,u'h,p'h^rh,Vh,qh) 

~ Bhi'^hi '^hy Phy Vh, Qh) + {{I Jh)^hy '^h)h dhTh)h 

“f {(^d JhjdhCTi^^ Vfi)ii “f ((/ Jfi^dhUf^^ dhViijh ((/ Jh)Phy ^h)h 

{{d Jh)Ui^, qh)h- 

Substituting the respective solutions (2.16) and (2.18) (and noting the slight discrepancy 
in the use of different harmonic forms), we have 

ddhi'^hy'^hy Phy'^hi Vfii Qh) — {'^hfy'^h)h {JhUfi,qh)h 

ddhi^^hy ^hy Phy ^hy ^hy Qh) {fhy^h)h {^hyQh^hy 

SO 

ddjiis^h ^hy '^hy Ph Phy '^hy '^hy Qh) 

= {fh - i*hfy ^h)h + «, qh)h - {why qh)h 

{,^d ((/ d}iTh)ii 

{(^d Jfi^djiCTV}i)ii ((/ Jh^dhUj^f dfiV}i)fi ((/ Jh^Phy ^h)h- 

As before, we bound the form above and below. For the upper bound, using Cauchy- 
Schwarz to estimate the extra inner product terms, we arrive at 

Bh{(Th -T,Uh- v,ph - q;Th,Vh,qh) 

< c (IIA - ilfWh + II - Wh)\\h + \\d - 411(11 Allv^ + IKIIvl + WPhWh) 
+\K - ^llu;, + IK - ^^Ili4 + \\p'h - QWh) (Ikhllu, + WvhWvn + \\qh\\h) ■ 

For the lower bound, we again choose (t/j, ah, qh) G to verify the inf-sup condition 
this time for Bh- 

Bh{(^h -T,Uh- v,ph - q;Th,Vh,qh) 

> Ih (I4h - r|4^ -f \\uh - ullu, + \\ph - q\\h) {\\rh\\v^ + \\vh\\v^ + ||gh|U) 

and 7 /i depends only on the Poincare constant cp||A|| ||7r/i||, uniformly bounded in h. 
Comparing with the upper bound and dividing out the common factor as before, this 
leads to: 

\Wh - r\\v^ -f \\uh - v\\v^ + \\ph - q\\h 

< C'jh^ iWfh - i*hf\\h + IIKA% - Wh)\\h +1|/ - 411(11 Allt7 + h'hWv^ + \\Ph\\h) 

+\Wh - '^Wv^ + IK - '>^\\v,, + \\p'h - (iWh) ■ 

Choosing (r, v, q) = (4, P^^p'y^), applying the triangle inequality with to account 

for the mismatch in the harmonic spaces, and using the well-posedness of the continuous 
problem (2.18), 

ll^/^ -<l|t4 + hh-u'y^Wv^ + \\Ph-p'h\\h 

< c (114 - 4/IU + WP^ini^h - Wh)\\h + 1|/ - 411(11/11 + KID + iIp^ - q\\h ). 

This differs from [22] in that we have the bound in terms of ||/|| + K||, and that we must 
estimate 114^^ (m)^ — Wh) ||h rather than ||45^“hlU alone. First, we use the modified Hodge 
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decomposition to uniquely write u'l^ as u'<g + Pfi'^u'h + with u'^ G and u'j^ G "5^^', 
and 

\\Pf^h{u'h-Wh)\\h < + M±)IU + 

(The projection is respect to the modified inner product). For the first term, we 
proceed exactly as in [22]: we have ~ ^ since the coboundary space is still the 

same, and thus only the term «'± contributes. Now u'j^ G so, using Jh to express it 
in terms of (/-orthogonality, we have Jhu'^ -L 3^, and thus Ps)f^Jhu'j_ = 0. Therefore, we 
have 

\\PM + «±)IU = \\P^,M\k = WP^.il - JhWJn <c\\l- AIKII/II + lliell). 

For the p'f^ term, this also proceeds as in [22] unchanged (except for, of course, the ex¬ 
tra 11^11 term): using the (unmodified) discrete Hodge decomposition, we have = 
P^^Ph + Pf-ihP'h = P^hP'h + *?• Since p’^ G S)'^, a similar argument gives Jhfi^ -L 53^, so 
P^JhP'h = 0 and 

lb; -q\\h = \\P<BAh = \\P<Bdi - AKIU <C\\I- MiWfW + IklD- 

Finally, we must consider the term \\Psj{Psj'^u'f^ — Wh)\\h- Expressing in terms of w, 
the terms do not combine as easily as the analogous terms involving fh and ilf, because 
their action as linear functionals operate on different harmonic spaces. 

Continuing with the proof of the theorem, we recall the third equation of (2.18): 

{Jhu'h,q')h = b>,b)h = {Jh{Jh%w),q')h 

which therefore says = P^/ i^w. This enables us to properly work with the 

modified orthogonal projection P^/^. Because is an isometry of the subspace bM4 to 
Wh , we have 

Ppjtw = itPi^^'W. 

where now is the orthogonal projection onto the identified image harmonic space. 
Then, using the triangle inequality again, 

WP^^APp^h -'u}h)\\h 

< ||p% -ilw)\\^ + \\Pf,fiJfi^ilw-ilw)\\h + \\Ps^S'i'hW-Wh)\\h 

< \\Psih\\ (INhIl II (^-^41-4)^11 + \\Jh^\\ 11 ^-All \\ilw\\h + \\i*hW - Wh\\h) 

< C (II (/ - P^i-yJ w|| + ||/ - All lAII + \\i*hW - Wh\\h) ■ 

The last term is the data approximation error for w, and the second term combines with 
the previous errors that reflect the non-unitarity of the operator. So, all that remains is to 
estimate the first term. Since it is in the subcomplex ihWh, the first estimate of Theorem 
2.11 applies: 

\\{I - PinpJ w\\ < \\{I - PMl < C inf llw-^lb, (2.24) 

by quasi-optimality. □ 

Concluding remarks of the proof of Theorem 2.16. To summarize, we proved our Main 
Theorem 2.16 by defining an intermediate solution on a modified complex that we iden¬ 
tify with a subcomplex, and analyzing the result via the Arnold, Falk, and Winther [3] 
framework. That theorem holds, with the estimate unchanged, though now u and Uh no 
longer are perpendicular to their respective harmonic spaces. The place where the extra 
terms all show up is in the variational crimes. In the process of arriving at a term that 
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looks like i^w — Wh, working with the different harmonic forms produces two more non- 
unitarity terms ||/ — Jh||(||/|| + HtfH), and finally, using Theorem 2.11 yields a direct 
estimate of how w fails to be a modified discrete harmonic form, giving the last best 
approximation term inf11 tu — ^ 11 y. □ 

We also note for future reference that in spaces where we have improved error esti¬ 
mates (which means tt/j are W -bounded maps) that we can replace that last V -norm in 
(2.24) to be the kk-inner product. Finally, we remark that, for a certain types of data 
interpolation, the errors \\fh — i*hf\\ and \\wh — ilMl be rewritten in terms of the other 
errors and another best approximation term. This will be useful for us in our examples. 

Theorem 2.18 (Holst and Stem [22], Theorem 3.12). IfUh ■ — >■ Ik^ is a family of 

linear projections uniformly bounded with respect to h, then for all f G kk^, 

||n,/ - 4/11 <c(\\I- 411 ll/ll + inf 11/ - /ll) . (2.25) 


3. Abstract Evolution Problems 

In order to solve and approximate linear evolution problems, we introduce the frame¬ 
work of Bochner spaces (also following Gillette, Holst, and Zhu [20]), which realizes 
time-dependent functions as curves in Banach spaces (which will correspond to spaces 
of spatially-dependent functions in our problem). We follow mostly [30] and [18] for 
this material. 


3.1. Overview of Bochner Spaces and Abstract Evolution Problems. Let X be a 

Banach space and / := [0, T] an interval in M with T > 0. We define 

C{I,X) := {u : I ^ X I u bounded and continuous}. 


In analogy to spaces of continuous, real-valued functions, we define a supremum norm 
on G(/, X), making C'(/, X) into a Banach space: 

\\u\\c{I,X) ■= sup ||M(f)||x- 

tGl 


We will of course need to deal with norms other than the supremum norm, which 
motivates us to define Bochner spaces: to define £^(/, X), we complete C{f X) with 
the norm 

r \ Vp 

\Ht)\\xdt] . 


\u\\lp{I,X) ■ — 




Similarly, we have the space X), the completion of C^{I, X) with the norm 

2 \ V2 


\U 


\\hhi,x)-.= / IKf)||^ + 


d , , 


dt 


X 


There are methods of formulating this in a more measure-theoretic way ([18, Appendix 
E]), considering Lebesgue-measurable subsets of /. 

As mentioned before, for our purposes, X will be some space of spatially-dependent 
functions, and the time-dependence is captured as being a curve in this function space 
(although this interpretation is only correct when we are considering C{I, X)—we must 
be careful about evaluating our functions at single points in time without an enclosing 
integral). Usually, X will be a space in some Hilbert complex, such as or 

where the forms are defined over a Riemannian manifold M. 
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We introduce this framework in order to be able to formulate parabolic problems more 
generally. It turns out to be useful to consider the concept of rigged Hilbert space or 
Gelfand triple, which consists of a triple of separable Banach spaces 

V CH CV* 

such that V is continuously and densely embedded in H. For example, if (V, d) is the 
domain complex of some Hilbert complex {W, d), setting V = and H = IF^ works, 
as well as various combinations of their products (so that we can use mixed formula¬ 
tions). H is also continuously embedded in V*. The standard isomorphism (given by the 
Riesz representation theorem) between V and V*, is not generally the composition of the 
inclusions, because the primary inner product of importance for weak formulations is the 
H-inner product. It coincides with the notion of distributions acting on test functions. 
Writing (•, •) for the inner product on H, the setup is designed so that when it happens 
that some F G V* is actually in H, we have F{v) = {F, v) (which is why we will of¬ 
ten write {F, v) to denote the action of F on r; even if F is not in H). In fact, in most 
cases of interest, the Tf-inner product is the restriction of a more general bilinear form 

between two spaces, in which elements of the left (acting) space are of less regularity 

than elements of H, while elements of the right space have more regularity. 

Given A G C{I, C{V, F*)), a time-dependent linear operator, we define the bilinear 
form 

a{t,u,v) ■.= {—A{t)u,v), (3.1) 

for {t,u,v) e M X F X F. To proceed, as in elliptic problems, we need a to satisfy some 
kind of coercivity condition, although it need not be as strong. It turns out that Garding’s 
Inequality is the right condition to use here: 

a(t,u,u) > Cilltilly — C 2 ||m||^, (3.2) 

with Cl, C 2 constants independent of f G /. Then the following problem is the abstract 
version of linear, parabolic problems: 

ut = A{t)u + f{t) 
m(0) = Uq. 

This problem is well-posed: 

Theorem 3.1 (Existence of Unique Solution to the Abstract Parabolic Problem, [30], 
Theorem 11.3). Let f G L‘^{I,V*) and uq G H, and a the time-dependent quadratic 
form in (3.1). Suppose (3.2) holds. Then the abstract parabolic problem (3.3) has a 
unique solution 

u G L\l,V)nH\l,V*). 

Moreover, the Sobolev embedding theorem implies u G C(I, H), which allows us to 
unambiguously evaluate the solution at time zero, so the initial condition makes sense, 
and the solution indeed satisfies it: m(0) = Uq. 

This theorem is proved via standard methods ([30, p. 382]); we take an orthonormal 
basis of H that is simultaneously orthogonal for F (a frequent situation occurring when, 
say, it is an orthonormal basis of eigenfunctions of the Laplace operator), formulate the 
problem in the finite-dimensional subspaces, and use a priori bounds on such solutions to 
extract a weakly convergent subsequence. With this framework, we can show that a wide 
class of PDE problems, particularly ones that are suited to finite element approximations, 
are well-posed. 


(3.3) 

(3.4) 
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3.2. Recasting the Problem as an Abstract Evolution Equation. Let us now see how 

these results apply in the ease of the Hodge heat equation (1.1) on manifolds. We take 
a slightly different approaeh from what is done in [20] and [1], solving an equivalent 
problem. This sets things up for our modified numerieal method detailed in later seetions. 

Let {W,d) be a elosed Hilbert complex, with domain complex {V,d), the standard 
setup in the above—in particular, we have the Poincare inequality and the well-posedness 
of the continuous Hodge Laplacian problem. We consider the space 2)^ •= x 1/^ 
and its dual x {V^)' with the obvious product norms (we use primes to 

denote dual spaces so as not to conflict with the dual complex with respect to the Hodge 
star defined earlier, though these uses are related). This, along with H = x W^, 

gives rigged Hilbert space structure 

The embeddings are dense and continuous by definition of the graph inner product and 
that the operators d have dense domain. We consider the Bocbner mixed weak para¬ 
bolic problem: to seek a weak solution {u,a) G L^(/,f5) ft to the mixed 

problem 

{a,uj) — {u,du) =0, f G/, 

{ut, (f) + {du, d(p) + {da, if) ={f,(p), tel, (3.5) 

m(0) = g, 

this makes it suitable for approximation using finite-dimensional subspaces of (e.g. 
degrees of freedom for finite element spaces). We see that (3.5) is the mixed form of 
(1.1), which amounts to defining a system of differential equations, introducing the vari¬ 
able cr defined by a = d*u, where d* is the adjoint of the operator d. We write the 
equation weakly (namely, moving d* back to the other side), which makes no difference 
at the continuous level, but will make a significant difference when discretizing. 

In order to use the abstract machinery above, we need a term with at. Formally differ¬ 
entiating the first equation of (1.2), and substituting = du in the second equation, we 
obtain 


0 = (atyOj) — (utydu) = {at,uj) — {f,du) + {da,du) + {du,ddu). 

Since d? = 0, that last term vanishes, and so, together with the equation for ut, we have 
the following system: 

{at,u!) + {da,du) = {f,du), tel, 

{ut, pi) + {da, if) + {du, dpi) ={f,p>), \/p>eV^, tel, (3.6) 

m(0) = g. 

Theorem 3.2. Suppose the initial condition g is in the domain of the adjoint V* and 
f G L^{I, {V^)'). Then the problem (3.6) is well-posed: there exists a unique solution 
{a,u) e L\I,^r))f^H\I,^')f^C{I,H) with (a(0),M(0)) = {d*g,g). 

Proof. We see that given / G Lf{I, {V^)'), we have that the functional F : {t,v) ^ 
{f,dT) -f {f,v) is in L'^{I,iT}'), since d maps to V^. For an initial condition on 
a, we can demand that (t( 0) be the unique ao statisfying (cto, r) — {g, dr) = 0. For this 

to reasonably hold, we must actually have at least uq G Vf, the domain of the adjoint 

operator d*, that is, ao = d*g. We equip the spaces with the standard inner products for 
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product spaces: 

{{a,u),{T,v))jj := (a,r) + {u,v) (3.7) 

(r,L)) 2 j := (a,r)y + {u,v)v (3.8) 

Consider the operator A :^ defined by 

a{a, u; u, ip) = {—A{a, u), (cc, ip)) = {da, du) + {da, ip) + {du, dtp). 

With the functional F defined as above, we have F G L‘^{1, 2)0’ so (3.6) is equivalent 
to the problem 

{a, u)t = A{a, u) + F. (3.9) 

We now need to verify that the bilinear form a satisfies Garding’s Inequality: 

a{a, u; a, u) = ||(icr||^ + {da, u) + \\du\\'^ 

— \W\\v ~ ll'^lP + {da,u) + llwlly — ||w||^ 

— Il'^lly “ ll'^ir ~ IM^II ll'^ll + ll'^lly ~ ll'^ll^ 

> Ikllv' - Ikll^ - \\W\\v - Okllv- + ll'“llv - Ikll^ 

= \ua,u)\\l-\\{a,u)\\P 

Thus, the abstract theory applies, and noting that the initial conditions {d*g, g) G H, we 
have that 

(a, u) G L\I, i^) n H\I, i^') n C{I, H) 

is the unique solution to (3.6) with initial conditions given by m(0) = g E and 
(j(0) = d*g. □ 

Given this, however, we must still establish that we also have a solution to the original 
mixed problem (which will be crucial in our error estimates): 

Theorem 3.3. Let {a, u) G LF‘{I, 2)) n H^{1, 2)0 n C{I, H) solve (3.6) above with the 
initial conditions. Then, in fact, (a,u) also solves (3.5). 

Proof. The second equation already holds, as it is incorporated unchanged into the equa¬ 
tions (3.6). To show the first equation, we show 

{at,u) - {ut,du) = 0 

for all time t. Then, since the original mixed equation holds at the initial time, standard 
uniqueness ensures it holds for all f G /. We simply realize it is setting p) = —dw. 

{at, u) - {ut, du) = {{a, u)t, {u, -du))H = a{at, up, u, -du) + (/, du) + (/, -du) 

= {da, du) + {da, —du) + {du, ddu) = 0. 

□ 


4. Error Estimates for the Abstract Parabolic Problem 

We now combine all the preceding abstract theory (the Hoist-Stem [22] framework 
recalled in §2.2, and the abstract evolution problems framework recalled in §3) to extend 
the error estimates of Gillette, Holst, and Zhu [20] and in particular, recover the case of 
approximating parabolic equations on compact, oriented' Riemannian hypersurfaces in 

'Using differential pseudoforms ([19, §2.8], [34]), we can eliminate this restriction. However, more 
theory needs to be developed for that case; the normal projection, in particular. We consider this in future 
work. 
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with triangulations in a tubular neighborhood. The key equation in the derivation 
of the estimates are the generalizations of Thomee’s evolution equations for the error 
terms. We shall see that these equations lead most naturally to the use of certain Bochner 
norms for the error estimates that are different for each component in the equation. 

Let {W,d) be a closed Hilbert complex with domain {V,d), and the Gelfand triple 
C Lf C on this complex as above. Now consider our previous standard setup 
of finite-dimensional approximating complexes {Wh, d) with domain (I 4 , d), with corre¬ 
sponding spaces = Vj^' X VI: (it is missing the harmonic factor), ih ■ Vh ^ V 
injective morphisms (that are VL-bounded), iih ■ Vh ^ V projection morphisms (which 
may be merely (/-bounded), and Wh o ih = id. Finally, we consider data interpolation 
operators II/i : VF —)■ Wh, such that Ilh o = id that realize which projections for the 
inhomogeneous and prescribed harmonic terms (fh and Wh in the abstract theory above) 
that we use. 

4.1. Discretization of the weak form. Suppose we have / G L^(/, iV^)') and g 
Let (<7,u) G L^(/,i5) n iF^(/,i5') n C{I,H) be the unique (continuous) solution to 
(3.5), as covered in §3. As in [20], we can consider approximations to this solution 
as functionals on finite-dimensional spaces e-g- finite element spaces. With the 
above considerations, we formulate the semi-discrete Bochner parabolic problem: 
Find {ah, Uh) I such that 


(UhiduJh^h 0 , 

{uh,t, ^h)h + {dah, (Ph)h + {duh, diph)h = {Vhf, (ph)h, 

Uh{0) = 9h- 


W Uh Vh 


tel 

tel 


(4.1) 

(We use the notation of Thomee for the test forms.) We define gh, the projected initial 
data, shortly. A similar argument as in §3 above, differentiating the first equation with 
respect to time, considering the Gelfand triple ^ I C Wj:-^ X W}: C (ig^)' gives that this 
problem is well-posed (or more simply, we choose bases and reduce to standard ODE 
theory as in (1.3) above). Following Gillette, Holst, and Zhu [20], we define the time- 
ignorant discrete problem, using the idea of elliptic projection [35] which we use to 
define a discrete solution via elliptic projection of the continuous solution at each time 
to e I: We seek {dh,Uh,Ph) £ ^h such that 

{^h, ^h)h (u/i, dUh} h 0, 

{ddh, (Ph)h + {duh, d(ph)h + {Ph, ^h)h = {Ilh{-Au{to)), (ph)h 

{uh,qh)h = {Ilh{Ps^u{to)),qh)h 


\/ Uh e 


k-l 


V G 


vi: 


y qh^ ^h- 


(4.2) 

Note that we have included a prescribed harmonic form given by the harmonic part of 
u (following [1]). We then take the initial data gh to be m/i( 0); it is just the solution to 
the elliptic problem with load data n/i(—A^f), since m(0) = g. Note we do not directly 
interpolate g itself via Yih for the data; the reason for this will be seen shortly. This 
discrete problem is well-posed, i.e., a unique solution Uh{to) always exists for every time 
to e /, by the first part of Theorem 2.16 above. The presence of an additional term 
ph and equation involving harmonic forms departs from Gillette, Holst, and Zhu [20], 
because the theory there is facilitated by the fact that there are no harmonic n-forms 
on open domains in MP (the natural boundary conditions for such spaces are Dirichlet 
boundary conditions, in contrast to the more classical example of 0-forms, i.e. functions). 
Here, however, we must consider harmonic forms, since we may not be working at the 
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end of an abstract Hilbert complex. For our model problem, namely differential forms 
on compact orientable manifolds (without boundary), even in the case of n-forms, the 
theory is completely symmetric (by Poincare duality [5, 25, 28]).^ In addition, the linear 
projections H/i may not preserve the harmonic space, which gives the possibility of a 
nonzero ph, despite — Am having zero harmonic part (so it is its own error term). 

4.2. Determining the error terms and their evolution. Continuing the method of 
Thomee [33], we use the time-ignorant discrete solution as an intermediate reference, 
and estimate the total errors by comparing to this reference and using the triangle in¬ 


equality. Roughly speaking, we try to estimate as follows: 

\\ih(^h{t) - cr(t)||y < \\ih(^h{t) - ih^h{t)\\v + \\ih^h{t) - cr(t)||y (4.3) 

\\ihUhit) - u{t)\\v < \\ihUh{t) - ihUhit)\\v + \\ihUhit) - u{t)\\v. (4.4) 

It turns out that this grouping of the terms is not the most natural for our purposes. We 
shall see it is the structure of the error evolution equations that groups the terms more 
naturally as: 

\\ihUhit) - u{t)\\ (4.5) 

\\ih(^h{t) - cr(t)|| -f \\d{ihUh{t) - u{t))\\ (4.6) 

\\d{ih(Th{t) - a{t))\\. (4.7) 


The sum of these three terms is the sum of the two (/-norms above. In addition, we shall 
see in our application to hypersurfaces that this particular grouping of the error terms also 
corresponds more precisely to the order of approximations in the improved estimates for 
the elliptic projection (namely, they are of orders and respectively, for 

degree-r polynomial differential forms). 

The plan is to use the theory of Holst and Stem [22] reviewed in §2.2 above to estimate 
the sum of the two second terms in (4.3) and (4.4); the elliptic projection simply is an 
approximation, at each fixed time, of the trivial case of u being the solution of the con¬ 
tinuous problem with data given by its own Laplacian, — Aw. The harmonic form portion 
will come up naturally as part of the calculuation. Using the notation of Thomee [33], 
we define the error functions 


p{t) := Uh{t) - ilu{t) 

(4.8) 

9{t) := Uhit) - Uh{t) 

(4.9) 

= ah{t) 

(4.10) 

e{t) := ah{t) - ah{t) 

(4.11) 


(Thomee does not define the third term -0; we have added it for convenience.) In the 
case that there are no variational crimes (i.e., Jh is unitary), the error terms p and xjj 
are bounded above by the elliptic projection errors (because there, is the orthogonal 
projection, and ||i^|| = ||i/i|| = 1), so that we have, for example, that \\ihUh — u\\ < 
||6*|| -f IIpII, corresponding to the use of p in [33, 20]. For our purposes, however, the 
choice of p here does not correspond as neatly, now being an intermediate quantity that 
helps us estimate 6 in terms the elliptic projection error (the second term in (4.4)). We 
find that it contributes more terms with \\I — Jh\\. Similar remarks apply for a and V’- We 
use the method of Thomee to estimate the terms 9 and £ in terms of (the time derivatives 

Despite this, there are a number of reasons why one should still prefer to continue to phrase problems 
in terms of n-forms if the problem calls for it ([19] describes how it affects the interpretation of certain 
quantities); and we shall see that it does in fact still make a difference at the discrete level. 



FEEC FOR EVOLUTION PROBLEMS ON MANIEOLDS 


23 


of) p and ijj, and the elliptic projection error; In order to do this, we need an analogue of 
Thomee’s error equations. 

Lemma 4.1 (Generalized Thomee error equations). Let 9, p, and e be defined as above. 
Then for all t E I, 

(e, Uh)h - {9, duh)h = 0 y(ph e 

{9t,(ph)h +{de,iph)h +{d9,d(ph)h = {-pt + Ph + (J^h - il)ut,(ph)h 'iuJhEV^. 

(4.12) 

This differs from Thomee [33] and Gillette, Holst, and Zhu [20] with the harmonic 
term ph, which accounts for the projections II/i possibly not sending the harmonic forms 
to the discrete harmonic forms, an extra d9 term which accounts for possibly working 
away from the end of the complex (for differential forms on an n-manifold, forms of 
degree k < n), and another data interpolation error term for Ut (which also distinguishes 
it from Arnold and Chen [1]). 

Proof. The first equation is simply weakly expressing e as df9. This follows immediately 
from the corresponding equations in the semidiscrete problem and the time-ignorant dis¬ 
crete problem. For the second term, consider the expression 

B := {9t,^h)h + {de,p>h)h + {d9,dph)h + {pt,Th)h, (4.13) 

and expand it using the definitions to obtain 

B Th)h Thjh 

-h {dah - ddh, (Ph)h + {duh - duh, dip)h + {uh,t, Th)h - {ilut, Th)h- 

We cancel the Uh,t terms, and apply the semidiscrete equation (4.1) to cancel the dah and 
duh terms, which gives us 

B = (rife/, iph)h - {ddh, Th)h - {duh, d(ph)h - Th)h, 

and finally, using the second equation of (4.2) to account for the middle terms, we have 

B = (rife/, ph)h + (nfe(AM), p)h + {ph, Th)h - {i*hUt, p>h)h 

= (rife (Am + f -Ut), (ph)h + {ph, Th)h + ((life - i*h)ut, (ph)h- 

But since Ut = Aw -|- / is the strong form of the equation, which we know is satisfied 
by the uniqueness, it follows that B = {ph + (life — iDut, ^h)h- Subtracting the pt from 
both sides gives the result. □ 

Now we present our main theorem. 

Theorem 4.2 (Main parabolic error estimates). Let {a, u) be the solution to the con¬ 
tinuous problem (3.5), {ah,Uh) be the semidiscrete solution (4.1), {dh,Uh) the elliptic 
projection (4.2), and the error quantities (4.8)-(4.11) be defined as above. Then we have 
the following error estimates: 

\\d{t)\\h < \\Pt\\L^I,Wh) + \\Ph\\L^I,Wh) + IKlIfe - i*h)Ut\\L^I,Wh) (4-14) 

\\d9(t)\\h + ||£(f)IU < C {\\pt\\L^(pWh) + \\Ph\\L^(I,Wh) + IKn^ - ih)'^t\\L‘2{I,Wh)) 

(4.15) 

(4.16) 


\\de(t)\\h < C [WftWL^^pWh) + IMh(II/i ~ Ph)'^t\\L‘^{i,Wh)) ) 
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with 


\\Pt\\L'^{I,Wh) — ^ {\\'i'hUh,t — Ut\\L'^{iy/) + ||/ — Jh\\c{Wh) (4.17) 

\\i’t\\L'^(I,Wh) ^ C {\\'ih^h,t — (^t\\L'^(I,W) + 11-^ — Jh\\c(Wh) \Wt\\L'^{ItW)) ■ (4.18) 

We may further eombine these terms, whieh we shall do in a separate eorollary, but it 
is useful to keep things separate, whieh allows terms to be analyzed individually when 
eonsidering speeifie ehoiees of V and I 4 . The error terms ihd'h — <J and ihiih — u and 
their time derivatives are furthermore estimated in terms of best approximation norms 
and variational erimes via the theory of Holst and Stern [22] . The different Boehner 
norms involved arise from the strueture of the error evolution equations. 

Proof. We adapt the proof teehnique in [33, 20] to our situation, and for ease of notation, 
unsubseripted norms will denote the PT-norms and norms subseripted with just h will 
denote norms on the approximating eomplex. We now assemble the estimates above 
separately by eomputing the PT-norms of the errors and their differentials. We begin by 
estimating ||6'(f) ||/j. We use the standard teehnique of using the solutions as their own test 
funetions: Set cph = 9 and = £ in (4.12). Adding the two equations together yields 

2+ Uplift + = (“Pi + + (n/i —(4.19) 

Following Thomee [33], we introduee 5 > 0 to account for non-differentiability at 6* = 0, 
and observe that 

= IjfPWl < (IlftlU + IIpoII* + ll(n<. - 

using (4.19), the Cauchy-Schwarz inequality, and the definition of operator norms (our 
goal is to get all of those quantities on the right side of the equation close to zero, so we 
need not care too much about their sign). Thus, since ||6'||/i <(l|9|II + i=)'A we have, 
canceling \\9\\h, 

f^m\l+< iiftiu + iimik + ii(n» - i’Mk- 

Now, using the Fundamental Theorem of Calculus, we integrate from 0 to f to get 

ll»(*)IU = ||»(0)||,. + limjf' J^ml + < /(IlftL + IIAII/. + IKn^- 

° ° (4.20) 

6*(0) vanishes by our choice of initial condition as the elliptic projection. 

Next, continuing to follow [20], we consider ||e:(f)||/j. We differentiate the first error 
equation and substitute iph = 29t and Uh = 2e, so that 

{et,2e)h - {6t,2de)h = 0 

{9t, 29t)h + {ds., 29t)h + {d9, 2d9t)h = {—pt + Ph + (n^ — iDut, 29t)h- 


(4.21) 

(4.22) 
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Adding the two equations as before, we have, by Cauehy-Sehwarz and the AM-GM 
inequality, 

j,\\4Um\\uyni 

— 2||pt ||h + 2||p/i||/i||0t ||/i + 2||(n/i — 

^ 2 (llpill^ + WPhWh + \\(J^h — ih)'^t\\h) + ill^tllL 

Again, dropping some positive terms (this time ||6*t||^), using the Fundamental Theorem 
of Calculus and noting the initial conditions vanish by the choice of elliptic projection, 
we have 

Iklli + ml < 2 f (lb,IIS + ||p»||S + ||(n„ - Il)u,\\l) . (4.23) 

Jo 

Finally, we estimate ||(i£||/i. As in the estimate above, we differentiate the first equation 
with respect to time, and substitute cu = 2et, = 2det, 


{et,2et)h- {0t,2det)h = 0 (4.24) 

{Ot, 2det)h + {de^ 2det)h + {dO^ 2ddet)h = {~pt + + (If/i — 2det)h- (4.25) 

Noting that = 0, is perpendicular to the coboundaries, and if) = d*f^p, we add the 
equations to get 

2\\£\\1 + = 2(-pi + (Rh - il)ut, det)h = 2{-'\\)t + dKUh - 4)«u £t)h 

< WP’tWh + IMKn/i - ih)'^t\\h + 2\\e\\l. 

By the Fundamental Theorem of Calculus, and noting vanishing initial conditions (and 
an exact cancellation of positive terms), we have 

ll<fe|lft< f‘ {\ml+\\dl{Tlk-r,)u,\\l). (4.26) 

Jo 

We now estimate p and ijj. We note that the time derivative of the solutions are also 
solutions to the mixed formulation, at least provided that Ut and other associated quan¬ 
tities are sufficiently regular (in the domain of the Laplace operator) for the norms and 
derivatives to make sense. Then (recalling we have 

\\p{t)\\h = \\uh - ilu\\ < \\uh - itu\\ + Wi^u - ilu\\ 

< IN/tll {\\ihUh - m|| + ||/ - Jh\\ ||m||) , (4.27) 
and 

\\i^{t)\\h = \\dh - ilaW < \\dh - iftCrll + ||f+cr - ila\\ 

< IKl II iW'i'hdh - ct|| + ll-f - Jh\\ Ikll). (4.28) 

The same estimates hold for the time derivatives. The first terms are the estimates that 
allow us to use the theory of §2.2. We note that the theory acutally uses L-norms, but 
it will work. We cannot improve this in the abstract theory; instead, we use theory for 
specific choices of V, W, and 14, such as appropriately chosen de Rham complexes and 
approximations to improve the estimates ([3, §3.5], [1, Theorem 3.1]). For these cases, 
it is helpful to keep the individual estimates on |4ip, etc. separated. We have 
combined terms because the abstract theory gives us all the variational crimes together, 
as it makes heavy use of the bilinear forms above. Additional improvement of estimates 
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based on regularity as done in [3] eannot made for the variational erimes, as diseussed in 
[22, §3.4]. We give the relevant example and result in the next seetion. □ 

Corollary 4.3 (Combined estimate). Let 9, p, ip, and e be as above. Then we have 

\\ihO'h — <T||Li(7,y) + \\ihUh — < 

c {\\pt\\L^{I,Wh) + ll(n/i - + \\'^t\\L^(pWh) + \\dh0^h - i*h)Ut\\L^{pWh) 

+ \\Ph\\L'^{i,Wh) + W'^hb'h — y) + \\ihUh — m||l2(/,v)) • (4.29) 
Further expanding the time derivative terms, we have 
— 0'\\L^(iy) + WihUh — uWii^y) < 

C {\\ihUh,t — y/) + \\ihd'h,t — (^t\\L'^{i,w) 

+ 11-^ — Jh\\ i|Mi||L2(/,m) + 11-^ — All ||(Tt||L2(/^vi/) 

+ ll(n/i - Ph)'^t\\L'^{I,Wh) + IIA(n/i - i*h)'^t\\L'^{I,Wh) 

+ \\ihPh\\L'^{I,W) + Wihdh — <T||L2(/^y) + \\ihUh — m||l2(/,V)) • 

These terms are organized as follows: the IC-error in the approximations of the time 
derivatives, the variational crimes with 11J—J/^ 11, the data approximation error for the time 
derivatives, and finally the V -approximation errors for the elliptic projection. These can 
be further expanded in terms of best approximation errors, but we will not have use for 
that outside of specific examples where the computation is easier done with the previous 
theorems. This corollary is simply stated for conceptual clarity and a qualitative sense of 
all the different individual contributions to the error. 

Proof. First, we note that by the Cauchy-Schwarz inequality, the estimate for Ijd^H (4.14) 
can be rewritten as using W) norms to match the squared terms (4.23) and (4.26). 
Combining and absorbing constants, we arrive at 

IAcT/i(f) — o'{t)\\v + WihUhit) — u{t)\\v < C {\\pt\\L'^(iy/h) + ll(nh — 'i*h)'^t\\L'^{i,Wh) 

+ \\ft\\L'^{I,Wh) + \\d-*h0^h - ih)'‘^t\\L^{I,Wh) + \\Ph\\L2{I,Wh)) 

+ lAA(f) - cr(f)||y + \\ihUh{t) - u{t)\\v. 

Integrating from 0 to T, the latter two C-norm terms become L^(I, V) norms (and ab¬ 
sorb the factor of T from integrating the first into the constant). Finally, using Cauchy- 
Schwarz to change the L^{I, V) norm into an L^(I, V) norm, and substituting for pt and 
ipt gives the result. □ 


5. Parabolic Equations on Compact Manifolds 

As an application of the preceding results, we return to our original motivating exam¬ 
ple of de Rham complex to explore an example with the Hodge heat equation on hyper¬ 
surfaces of Euclidean space, generalizing the discussion in [22, 20]. Let M be compact 
hypersurface embedded in M inherits a Riemannian metric from the Euclidean 

metric of 


5.1. The de Rham Complex on a Manifold. We define the Lf differential fc-forms on 
M given by 



ay.ydx"^ A • • • A G 


Hi...Ik 


G A(M) 
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the standard indexing of differential form basis elements, namely strictly increasing se¬ 
quences from {1,..., n}. The inner product is given by {u,r]) = J to A *17, where a is 
the Hodge operator corresponding to the metric. 

The weak exterior derivative is defined on the domains HVL^{M), and we have a 
Hilbert complex d) with domain complex d), with d^^^ o d^ = 0: 

0-^ ^ 0. 

As required in the abstract Hilbert complex theory, each domain space carries the graph 
inner product: 

{u,v )•= {u,v)L2^ki^M) + {d'^U, d^v)L2Qk+i(^]^y 

For open subsets U C the ends (fc = 0 and k = n) of this complex are familiar 
Sobolev spaces of vector fields with the traditional gradient, curl, and divergence opera¬ 
tors of vector analysis: 

0-- H\U) H{cnT\) -- i7(div) L\U) -- 0. 

Similarly, the dual complex is H*fl{M) defined by := con¬ 

sisting of Hodge duals of (n — fc) -forms. We have that the embedding HVL^{M) ft 
H*VL^{M) C‘^VL^{M) is compact, which enables a Poincare Inequality to hold and 
the resulting Hilbert complex d) to be a closed complex [29, 3]. To summa¬ 

rize, we have the following: 

Theorem 5.1. Let M be a compact Riemannian hypersurface in Then taking 

with maps df the exterior derivative defined on the domains = 
HQfi{M), (kF, d) is a closed Hilbert complex with domain (F, d). 

We thus are able to define Hodge Laplacians, and see all the abstract theory for the 
continuous problems (2.15) and (3.6) applies with these choices of spaces. 

5.2. Approximation of a hypersurface in a tubular neighborhood. In order to ap¬ 
proximate the problems (2.15) and (3.6), we consider, following [22], a family of ap¬ 
proximating hypersurfaces to an oriented hypersurface M all contained in a tubular 
neighborhood U of M. The surfaces Mh generally will be piecewise polynomial (say, 
of degree s); the case s = 1 corresponds to (piecewise linear) triangulations, studied in 
[16, 14], and generalized for s > 1 in [13]. However, the piecewise linear case still is 
instrumental in the analysis and indeed, the definition of the spaces (via Lagrange inter¬ 
polation), and so we shall denote it by Th (the triangulation, i.e., set of simplices, will be 
correspondingly denoted by Th, and their images under the interpolation will be denoted 
Th)- It is convenient, also, to assume that the vertices of the both the triangulation and 
the higher-degree interpolated surfaces actually lie on the true hypersurface. 

The normal vector u to the M allows us to define a signed distance function dist : 
U M. given by 

dist(x) = ± dist(x, M) = ± inf \x — y\ 

y£M 

where the sign is chosen in accordance to which side of the normal x lies on. By ele¬ 
mentary theorems in Riemannian geometry [15, Ch. 6], dist is smooth, provided U is 
small enough; the maximum distance for which it exists is controlled by the sectional 
curvature of M. The normal a can be extended to the whole neighborhood; in fact it is 
the gradient V(5. It is also convenient to define the normals Uh to the approximating sur¬ 
faces Mh- In most of the examples we consider, we assume the vertices of Mh (and Th) 
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U 



Figure 1. A curve M with a triangulation (blue polygonal curve M^) 
within a tubular neighborhood U of M. Some normal vectors u are drawn, 
in red; the distance function 5 is measured along this normal. The inter¬ 
section X of the normal with defines a mapping a from x to its base 
point a{x) G M. 

lie on M, but this is not a strict requirement. Instead, we need a condition to ensure that 
the hypersurfaces Mh are diffeomorphic to M, eliminating the possibility of a double 
covering (e.g., as pictured in [17, Fig. 1, p. 12]). In particular, we want Mh to have the 
same topology as M. This is again restriction on the size of the tubular neighborhood. 
In such a neighborhood U, every x e U decomposes uniquely as 

a; = a(a:) + dist(a:)i/(a:), (5.1) 

where a{x) G M, and a : f/ —)■ M is in fact a smooth function, called the normal pro¬ 
jection. a can then be used to define the degree-s Lagrange interpolated hypersurfaces 
by considering the image of Th under the degree-s Lagrange interpolation of a over each 
simplex in Th (we write Uk '■ Th ^ Mh for this) [13, §2.3]. Now, Holst and Stem [22] 
show, for hypersurfaces, the following result for the variational crime || J — J/i||: 

Theorem 5.2 (Holst and Stern [22], Theorem 4.4). Let M be an oriented, compact 
m-dimensional hypersurface in and Mh be a family of hypersurfaces lying in 

a tubular neighborhood U of M transverse to its fibers, such that ||()||oo —^ 0 and 
— T^hWoo — >■ 0 as /i —)■ 0. Then for sufficiently small h, 

\\I-Jh\\<Cm\oo + \W-i^h\\l). (5.2) 

A result of Demlow [13, Proposition 2.3] states that, in the case that Mh is obtained 
by degree-s Lagrange interpolation, that ||(5||oo < and \\u — UhWoo < Ch^. Thus, 

putting these results together, we have that 

\\I-Jh\\<Ch^+\ (5.3) 

Now, the three best approximation error terms (2.17) for finite element approximation 
by polynomials of degree r are bounded by CT, or CT~^, depending on the 

component chosen, so it is crucial to allow for this case, and the convergence rate is 
optimal when r = s. Figure 2 also dramatically demonstrates how much better a higher- 
order approximation can be with a given mesh size. 
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Figure 2. Approximation of a quarter unit circle (black) with a segment 
(blue) and much better quadratic Lagrange interpolation for the normal 
projection (red), though the underlying triangulation is the same (and thus 
also the mesh size). 


Restricting a to the surfaces M^, gives diffeomorphisms 

a : Mh M is therefore a diffeomorphism when restricted to each polyhedron (and is at 
least globally Lipschitz continuous, the maximum degree of regularity in the piecewise 
linear case. This is not a problem for Hodge theory, because the form spaces are at most 
where regularity is concerned; see [36]). See Figure 1. 

5.3. Finite element spaces. We choose finite-dimensional subspaces of HVL^{Mh) 
for each k, satisfying the subcomplex property d/iA^ C We can then pull forms 

on Mh back to forms on M via the inverse of the normal projection, which furnishes the 
injective morphisms : A^ HVL^{M) (since pullbacks commute with d) required 
by the theory above in Section 2. The main finite element spaces relevant for our pur¬ 
poses are two families of piecewise polynomials, discussed in detail in [2, 3]. We must 
choose these spaces for our equations in a specific relationship in order for the numerical 
methods and theory detailed above to apply, and for the approximations to work. This 
is why we prefer a piecewise polynomial approximation of M as opposed to a curved 
triangulation of M itself; these are shown to have these necessary properties. 

Definition 5.3 (Polynomial differential forms). Let Vr denote polynomials of degree at 
most r, in n variables, and LLr be the subspace of homogeneous polynomials. We define 
the first family, denoted VrA.^{T), to consist of all k-forms with coefficients belonging 
to Vr when restricted to each n-simplex of T. The continuity condition is that the poly¬ 
nomials on two simplices having a common face must have the same trace to that face. 
The second family, denoted VfA^{T), are intermediate spaces, between the spaces of 
the first class: 


Vr-iA\T) C V;K\T) C VrA\T). 
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These are defined as follows: first, consider the radial vector field X = that is, at 

each X, it is a radially pointing vector of length |a;|, and then define the Koszul operator 
Ku \= X A u, the interior product with X. Then 

V;A\T) := 

This is a direct sum, since k always raises polynomial degree and decreases form degree, 
so yields homogeneous polynomials of degree r. k is in some ways dual to the operator 
d (which, in particular, increases form degree and decreases polynomial degree), and by 
the properties of interior products, rf = 0. 

These polynomial spaces generalize existing finite element spaces, such as Whitney 
forms, Nedelec elements, and Raviart-Thomas elements (see [20, 3] for these examples 
and more), realizing the collection and clarification of previous results respecting vector 
methods, as we have mentioned numerous times throughout this work. 

The important property of these spaces is that they admit the cochain projections 
whose role we have seen is so important in the theory. First, we describe the case where 
M = U is a domain in MT with smooth or Lipschitz boundary. 

4 : kl where e {Vrk^(T), V;A.^(T)}. (5.4) 

These operators, by virtue of their construction, are uniformly bounded (in LfOf, not 
just HOf) with respect to h. Finally, the following theorem explicitly expresses the 
projection error (and hence, best approximation error) in terms of powers of the mesh 
size h and the norms of the solution. 

Theorem 5.4 (Arnold, Falk, and Winther [3], Theorem 5.9). 

(i.) Let be one of the spaces V~^_iA^(T) or, if r > 1, Vrk^iT). Then vr^ is a 
cochain projection onto A^ and satisfies 

for 0 < s < r + 1. Moreover, for all u G LfOffJ), tt^o; —)■ u in L'^ as h ^ 0. 
(ii.) Let A’l be one of the spaces VrA^(T) or Vf A\T) with r > 1 . Then 

\\d{u — Tih^)\\L2nt‘{u) < ch^\\duj\\H‘n’^{u)^ ^ ^ iF^f 2 ^(f/), 

for 0 < s < r. 

These bounded cochain operators are explicitly constructed in [2, 3]; they are the 
natural interpolation operators defined for continuous differential forms and analogous 
to polynomial interpolation operators on functions, but combined with smoothings to 
allow extension to differential forms which may not necessarily be continuous. 

Example 5.5 (The Mixed Hodge Laplacian problem on an open subset of M"^). For the 
mixed Hodge Laplacian problem we considered above, we must choose A^“^ and A’j in 
such a manner such that dA’jf^ C A|; one cannot make the choices of spaces completely 
independent of one another for our mixed problem [3, §5.2]. For example, if we choose 
Al-^ = VrA^-\Th), we necessarily must choose 

Ale{Vr-iA\rH), V;A\rH)]. 

Similarly, for A^~^ = Vf A^~^(Th), we choose 

Ale{V;A\TH), Vr-iA\Th)]. 

Continuing in this manner down the complex, there are 2” possible full cochain sub¬ 
complexes one can form with these choices of spaces. Of course, for one single Hodge 
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Laplacian problem, we only need to work with three spaces in the chain, since the equa¬ 
tions only involve {k — 1)- and k-forms and their differentials. 

Example 5.6 (Finite Element Spaces on Riemannian manifolds). Now, suppose we are 
back in the situation with a Riemannian hypersurface M C with a family of 

degrees Lagrange-interpolated surfaces M^, over a triangulation T^. We can still con¬ 
sider the polynomial finite element spaces on the triangulation Th as before; the only 
difference here is that the simplices may not join up smoothly (i.e., as a manifold, it may 
have corners). This is not a problem, because the continuity conditions enforced by the 
finite element spaces also allow for discontinuities or non-classical-differentiability on 
the simplicial boundary faces. To define the analogous polynomial spaces on the pos¬ 
sibly curved triangulations M^, we simply say a form is in the analogous polynomial 
spaces VrN^iTh) if its pullback by the inverse of the interpolated normal projections 
Ok '.Th ^ Mh to Th is in VrA^{Th) [13, §2.5]. Now, from VrJk^iTh), we pull these forms 
back to the surface M via the normal projections This gives the injective 

morphisms : A’j ^ HQ^{M); it commutes with the differentials, since the pullbacks 
do. 

For the bounded cochain operators, the situation is similar. We haven'll : HQfi[Mh) —)■ 
A.\ a cochain projection defined by pulling forms defined in neighborhoods back to the 
triangulations (using the trace theorem if necessary), as constructed in [2, 3]. Then 
we compose with the pullbacks This gives us the cochain projections : 

HQfi(M) —)■ A| (by [22, Theorem 3.7]). 

5.4. Estimates for the Mixed Hodge Laplacian problem on manifolds. With this, we 
can then integrate the terms from [22, Example 4.6] to get the results for the parabolic 
equations (or, equivalently, add the variational crimes to [20, 1]). Let us consider now the 
mixed Hodge Laplacian problem on Riemannian hypersurfaces, considering the setup in 
the previous example. Namely, we consider = HQf(M) as above, 

the approximating spaces = Vr+iA^~^(Th) and = VrN^{Th), and finally the 

inclusion and projection morphisms as above (possibly with additional pullbacks for 
interpolation degree s > 1). Of course, as mentioned before, these are not the only 
ways of choosing the spaces, but we stay with, and make estimates based on, this choice 
for the remainder of this example (the same choice made in [22, Example 4.6]). For a 
function / G LF‘Qf{M), we have an approximate solution pj) G ih'^'h ^ the 

elliptic problem, on the true subcomplex ihWh (with modified inner product, as in the 
theory of §2.4). For / sufficiently regular, and (a, u,p) satisfying the regularity estimate 
[3, 20] 

||w||//s+2 + ||p||jy^+2 + ||cr||H^+l < C||/||//S, (5.5) 

for 0 < s < Smax, then, since we are in the de Rham complex, where the cochain 
projections are VF-bounded, we have the improved error estimates of Arnold, Falk, and 
Winther [3, §3.5 and p. 342] for the elliptic problem: 

11“ - ihu'hW + lb - ^hP'hW < CN^^WfWnr-i (5.6) 

\\d(u - ihUh)\\ + ||cr - iha'hW < CN\\f\\Hr-i (5.7) 

\\dia-tha',)\\<CN-^f\\Hr-i. (5.8) 

We should also note that Arnold and Chen [1] prove that this also works for a nonzero 
harmonic part [1, Theorem 3.1]. Holst and Stem [22] augment these estimates to in¬ 
clude the variational crimes, so that (changing the notation to suit our problem) for 
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(cT/i, Uh, Ph) e Xh, the discrete solution to the elliptic problem now on the approximating 
complexes we have chosen, we have the estimates 

IL - ihiihW + Up - ihPhW + - ih&h)\\ + Ik - ikikW) 

+ - ikffkil < C(A’'+‘||/1L-. + A*+'ll/ll). (5.9) 

We note the terms associated to the different powers of h above correspond exactly to 
the breakdown (4.14)-(4.16) above. For the elliptic projection in our problem, we also 
need to account for the nonzero harmonic part of the solution. Setting w = Ps^u and 
Wh = If/itf}, we have that our three additional terms (given by Theorem 2.16 above) are 
the corresponding best approximation error ini^^yk ||ti} — t>||v, the ||/ — J/i|| term, and 
the data approximation \\wh — ihw\\h- For the best approximation, we make use of our 
observation about the inequality (2.24), in which we may instead use the PF-norm instead 
of the F-norm in the case that the projections are IF-bounded, as they are here in the de 
Rham complex. Because w is harmonic, it is smooth (and in particular, in so we 

may apply Theorem 5.4 to find that it is of order C 11 te 11 h^+i . The 11 / — A11 term has 
already been shown to be of order above in Theorem 5.2. Finally, by Theorem 

2.18 above, we have that data approximation splits into the other two terms. Therefore, 
to summarize, we have 

Theorem 5.7 (Estimates for the elliptic projection). Consider the solution 

to the parabolic problem (3.6) and {ahit),Uhit)) the semidiscrete solution in (4.1) above. 

Then we have the following estimates for the elliptic projection (d/i, Uh,ph)-' 

\\u - ihUhW + \\ihPh\\ + h {\\d{u - ihUh)\\ + lk - ih^hW) 

+ h‘^\\d{a - ihdh)\\ < C (||AM||j^r-i + ||w}||iyr+i) + (||Am|| + \\w\\)) . 

(5.10) 

(We note p = P^f—Au) = 0.) We now would like use the our main parabolic esti¬ 
mates to analyze the analogous quantity 

\\u{t)-ihUh{t)\\+h{\\d{u{t) -ihUh{t))\\ + \\a{t) - ih(Th{t)\\) + h‘^\\d{a{t) -ihah{t))\\, 

(5.11) 

and its integral, i.e. Bochner norm. 

Theorem 5.8 (Main combined error estimates for Riemannian hypersurfaces). Let {cr{t ), u{t)), 
{ahit),Uhit)), and all terms involving the elliptic projection are defined as above, and 
the regularity estimate (5.5) is satisfied. Then 

||w — ihUh\\L^(W) 

+ h (||ci(M - ihUh)\\L^{w) + Ik - AO-hlUqw)) + h^\\d{(T - ih(yh)\\L^{w) 

< C ((T + 1) (II AM||il(//r-l) -f ||'Ui||Ll(//r+l)) + T (II Arttll^qj^r-l) -I- ||'Uii||Ll(//r+l))) 

+ ((T + 1) (IIAm||2,i(w) + IkllLi(w)) + T (IIAm* 11^1 (VI/) -f IlfiijLqw)))] • 

(We abbreviate X) as L^{X).) The constants T, of course, can be further rolled 
into the constant C. We remark that in previous results, factors of T show up on the 
II AmJ terms, and, heuristically speaking, this is due to the Ut being a physically different 
quantity, namely, a rate of change. However, the appearance of the factor of T on the 
II Am|| comes from the harmonic approximation error which is, physically speaking, a 
harmonic source term. The details depend on the nature of the approximation operators 
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Proof. By the triangle inequality, we have that (5.11) breaks up into something of the 
form (5.9) (taking {dh, Uh, Ph) to be elliptic projection with / = —Au{t) and p = 0; here 
/ is not to be confused with the parabolic source term /(f)) and 

||4|| {\\9{t)\\h + h{\\e{t)\\h + \\du{t)\\h) + h^\\de(t)\\h) , (5.12) 

recalling the error quantities defined in (4.8)-(4.11). Now, substituting our estimates 
(4.14)-(4.16), we then have 

\\9(t)\\h < llptWi^Wh) + WPhWLHWh) + ll(n/i - ih)'^t\\LHWh) 

< C {\\ihUh,t — UtWi^iw) + IIp/iIUmw;,) + ||/ — Jh\\ IlMtlUqw) + ||(nh — ih)'^t\\L^{Wh)) 

< (II+ \\Aut\\L^{H^-^) + ||w||l1(_H'’'+1) + ||Wt||Li(iT’'+i)) 

+ C2h^~^^ [\\Au\\l^(^w) + ||^'*^t||Li(w) + Il'M^llLqw) + Il'iZ’tlUqw)) • (5.13) 

For \\d9\\h + II^IU, the computation is almost exactly the same, except with possibly 
different constants, to account for using Bochner norms, and that: 

< C [\\ihUh,t — Ut\\L-^(w) + IIPl||l2(W;,) + 11-^ “ Jh\\ ||Mt||L2(^y) + 11(11/1 — Ph)'^t\\L'^{Wh)) 

(II AM||i 2 (J^r-l) + II AMt||i 2 (J^r-l) + 11 fi) 11 L 2 (J^r +1 ) + 11 W)/11 ^2 r+ 1 ) ) 

+ (II AM||2,2(^y) + II AMi||2.2(W) + 11 11^2(14/) + ||'WijL2('4^/)) . 

These terms are actually absorbed into the lower order terms by the extra factor of h, due 
to consisting entirely of the same order terms except using a different norm. However, 
the situation is slightly different for ||(i£||/i; namely we use (5.7) to get a term of order h'", 
and the d^ on the variational crime part also removing a factor of h: 

\\de{t)\\h < C {\\ft\\L‘ 2 (Wh) + \\d*h(n.h - 'ih)'^t\\L^iWh)) 

< C {\\ihb'h,t — 0't\\L2(w) + 11-^ — Jh\\ ||(Tt||L2(44/) + \\dl_{Ilh — ih)'’^t\\L‘2{Wh)) 

< C^h'' {\\Aut\\L'2(^H'"-i) + ||M)||£, 2 (J/-r+l)) 

+ C'6^^(||^'w||L2(vy) + II Am/11^2 (44/) + ||m)||l 2(W) + ll'^^'t ||L2(vy)) 

However, we see that multiplying by hf, this term also gets absorbed; thus we need only 
consider the error from Hd^H/i in further calculation of the combined estimate. We have, 
thus far: 

||M(f) - ihUh{t)\\ 

+ h (||(i(M(f) - ihUhit))\\ + ||(T(f) - ihO-hit)\\) + li^Wdiait) - 4cr/j(f))|| 

< (IIAM||il(//-r-l) + IIAMt||j;,l(//-r-l) + ||w)||Li(_H-r-+i) + ll'lZ'tllLqHr'+i)) 

+ C2h^~^^{\\Au\\Li(^w) + II^'wJlI(W) + II'?^I|l1(W) + II'*^Jl1 (W)) 

+ C (||AK(f)||^/r-l + \\w{t)\\Hr+l) + + ^ (||AM(f)|| + ||w(f)||)) . (5.14) 

Integrating with respect to t from 0 to T, we find that the already-present Bochner norms 
are constant and thus introduce an extra factor of T. Absorbing the constants except T 
gives the result. □ 

This shows, in particular, that the optimal rate of convergence occurs when r = s, i.e., 
the polynomial degree of the finite element functions matches the degree of polynomials 
used to approximate the hypersurface. This tells us, for example, it is not beneficial to 
use higher-order finite elements on, say, a piecewise linear triangulation. Finally, to put 
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these estimates into some perspeetive and help develop some intuition for their meaning, 
we present the generalization of the estimates of Thomee from the introduetion. 

Corollary 5.9 (Generalization of [33, 20, 1]). Focusing on just the components u and 
a separately, we have the following estimates (assuming the regularity estimates (5.5) 
are satisfied), and supposing r = s, i.e., the finite element spaces considered consist of 
polynomials of the same degree as the interpolation on the surface: 

||M(t) -ihUh{t)\\ < (\\u{t)\\Hr+i + [ (||M(s)||jyr+i + ||Mt(s)||^r+i) ds 




(t) -ih(^h{t)\\ < M|M(f)||j^r +2 + (||M(s)||^r+l + \\Ut{s)\\Hr+l) ds^ 


1/2N 


This easily leads to an estimate in a Boehner L°° norm (simply take the sup in the 
non-Boehner norm terms and t = T in the integrals); this shows that the error in time is 
small at every t G I. Similar estimates hold for L^(/, W) norms. 

Proof We eonsider the improved error estimate and variational erimes in u and a sepa¬ 
rately. We first have, by expanding the terms in (4.20) as in the derivation of (5.13), 

\\u{t) - ihUhit)\\ < C{\\u{t) - ihUh{t)\\ + ||6'(f)||) 

< (\\Au{t)\\Hr-i + ||w(t)||/fr+l 


+ / (||AM(s)||//r-i -f ||AMt(s)||H-i + ||te(s)||H-i + \\u!tis) 




ds 


The result follows by noting that ||M||j|^r+i ineludes estimates on all the seeond derivative 
terms in u, and w = P^u, so those two norms ean all be eombined (with possibly differ¬ 
ent eonstants). Next, we eonsider a. The improved error estimates [3, p. 342] imply that 
if we do not combine estimates involving du with those of a for the modified solution, 
and / is regular enough to use the H^- rather than iJ^“^-norm, then we can gain back 
one factor of h, so that it is of order (rather than h'" as in (5.7)). On the other hand, 
the elliptic projection error ||£(f) || still can be taken along with \\da{t) || and was of order 
to begin with. Thus, applying (4.23), we have 

||(T(f) - ihcrh{t)\\ < C (||cr(f) - ihdh{t)\\ + ||£(f)|| + |MM(f)||) 


<Ch^+^ ||AM(f)||H^ + ||fi} 


ffr+l 


Uo 


(II AM(s)||^r-l -f II AMt(s) ll^r-l + ||'?T'(s) ll^r-1 + 11 (s) 11 H’’” ^ 


< ( \\u{t)\\Hr +2 + ( / (||m(s)||^.+i -f ||Mt(s)||^r+i) ds 


1/2N 


1/2N 


where we have used the same consolidation techniques for the norms on Au and w into 
norms on u as before. □ 


We see the variational crimes (arising from the extra pf) account for the sole additional 
term in the integrals. This cannot be improved without further information on the pro¬ 
jections rift. Otherwise, for r = 1, which correspond to piecewise linear discontinuous 
elements for 2 -forms (u), and piecewise quadratic elements for 1 -forms (a) with normal 




FEEC FOR EVOLUTION PROBLEMS ON MANIEOLDS 


35 


continuity (Raviart-Thomas elements), as studied by Thomee, we obtain the estimates he 
derived (and since the ph is not there in his case, we have that the extra terms with u do 
not appear under the integral sign). 


6. A Numerical Example 


In order to actually simulate a solution to the Hodge heat equation, we consider the 
scalar heat equation on a domain in M C but now using a mixed method with 2 - 
forms rather than the functions. We return to the evolution equation for both a and u, 
(3.6) above, which we recall here: 


tel, 

tel. 


( 6 . 1 ) 


{at,uj) + {da,du) ={f,du), \/uj eV’^ 

{ut,^p) + {da,ip) + {du,dip) ={f,^), 'iipeV’^, 

m(0 ) = g. 

Given Sh = HVL^[M) and = HVP{M), we choose bases, and use 

the semidiscrete equations (1.4), which we recall here (setting U to be the coefficients of 
Uh in the basis for Sh, and S to be the coefficients of ah in the basis for Hh) 

d fD -B'^\ / 0 0 \ /o' 


dt \0 A 


U 


-B -K 


U 


+ 


( 6 . 2 ) 


This may be discretized via standard methods for ODEs. Eor our implementation, we 
use the backward Euler method. This means we consider sequences (S"^, U"') in time, 
and then rewrite the derivative instead as a finite difference, evaluating the vector field 

portion on the right side at timestep n + 1, taking M = 



or 

;))(,".“)=-“(y-“(A) 

We now have written the system as a sparse matrix times the unknown f/”+^). 

This allows us to solve the system directly using sparse matrix algorithms without ex¬ 
plicitly inverting any matrices, making the iterations efficient. To analyze the error of 
the approximations, we can combine the above error estimates with the standard error 
analysis of Euler methods. See Eigure 3. 


7. Conclusion and Euture Directions 

We have seen that the abstract theory of Hilbert complexes, as detailed by Arnold, 
Ealk, and Winther [3], and Bochner spaces, as detailed in Gillette and Holst [20] and 
Arnold and Chen [1], has been very useful in clarifying the important aspects of elliptic 
and parabolic equations. The mixed formulation gives great insight into questions of 
existence, uniqueness, and stability of the numerical methods (linked by the cochain 
projections tt/j). The method of Thomee [33] allows us to leverage the existing theory 
for elliptic problems to apply to parabolic problems, taking care of the remaining error 
terms by the use of differential inequalities and Gronwall estimates (in the important 
error evolution equations (4.12) above). Incorporating the analysis of variational crimes 
allow us to carry this theory over to the case of surfaces and their approximations. 

We remark on some possible future directions for this work. Some existing surface 
finite elements for parabolic equations have been studied by Dziuk and Elliott [17], and 
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(C) 3 seconds (d) 4 seconds 


Figure 3. Hodge heat equation for A; = 2 in a square modeled as a 100 x 
100 mesh, using the mixed finite element method given above. Initial data 
is given as the (diseontinuous) eharaeteristie funetion of a C-shaped set in 
the square. The timestepping method is given by the baekward Euler 
diseretization, with timestep At = 5 x 10“®. 


mueh other work by Dziuk, Elliott, Deekelniek [12, 11], whieh aetually treat the case 
of an evolving surface, and treat a nonlinear equation, the mean curvature flow. Gener¬ 
ally speaking, this translates to an additional time dependence for evolving metric coeffi¬ 
cients, and a logical place to start is in the Thomee error evolution equations (4.12). Non¬ 
linear evolution equations for evolving metrics also suggests the Ricci flow [27, 9, 10], 
instrumental in showing the Poincare conjecture. The challenge there, besides nonlin¬ 
earity, is that tensor equations do not necessarily fit in the framework for EEEC. On the 
other hand, the Yamabe flow [31], which solves for a conformal factor for the metric (and 
is equivalent to the Ricci flow in dimension 2) suggests an interesting nonlinear scalar 
evolution equation for which this analysis may be useful. 

Gillette, Holst, and Zhu [20] also analyzed hyperbolic equations in this framework, 
and it would be interesting and useful to analyze methods on surfaces (including the 
evolving case), as well as taking a more integrated approach in spacetime. This is usually 
taken care of using the discrete exterior calculus (DEC), the finite-difference counterpart 
to EEEC to analyze hyperbolic equations [26] . 
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